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^ \ We have analyzed FUSE (905-1187 A) spectra of a sample of 16 local star- 

burst galaxies. These galaxies cover almost three orders of magnitude in star- 
CN \ formation rates and over two orders of magnitude in stellar mass. Absorption 

features from the stars and interstellar medium are observed in all the spec- 
tra. The strongest interstellar absorption features are generally blue-shifted by 
~ 50 to 300 km s~^, implying the almost ubiquitous presence of starburst-driven 
galactic winds in this sample. The outflow velocites increase with both the star 
^ ' formation rate and the star formation rate per unit stellar mass, consistent with 

a galactic wind driven by the population of massive stars. We find outflowing 
^ coronal-phase gas (T ~ 10^'^ K) detected via the O VI absorption-line in nearly 

00 ' every galaxy. The O VI absorption-line profile is optically-thin, is generally weak 

^ near the galaxy systemic velocity, and has a higher mean outflow velocity than 

' seen in the lower ionization lines. The relationship between the line width and 

column density for the O VI absorbing gas is in good agreement with expecta- 

■ tions for radiatively cooling and outflowing gas. Such gas will be created in the 
O ■ interaction of the hot out-rushing wind seen in X-ray emission and cool dense 

ambient material. O VI emission is not generally detected in our sample, sug- 
^ ■ gesting that radiative cooling by the coronal gas is not dynamically significant in 

■ draining energy from galactic winds. We find that the measured outflow veloc- 
ities in the HI and HII phases of the interstellar gas in a given galaxy increase 
with the strength (equivalent width) of the absorption feature and not with the 
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ionization potential of the species. The strong hnes often have profiles consist- 
ing of a broad and optically-thick component centered near the galaxy systemic 
velocity and weaker but highly blue-shifted absorption. This suggests that the 
outflowing gas with high velocity has a lower column density than the more qui- 
escent gas, and can only be readily detected in the strongest absorption lines. 
From direct observations below the Lyman edge in the galaxy rest-frame, we find 
no evidence of Lyman continuum radiation escaping from any of the galaxies 
in the sample. Moreover, the small escape fraction of light in the center of the 
strong C II absorption feature confirms the high opacity below the Lyman limit 
in the neutral ISM. The absolute fraction of escaping Lyman continuum photons 
is typically <1%. This sample provides a unique window on the global proper- 
ties of local star-forming galaxies as observed in the far-UV and also provides a 
useful comparison sample for understanding spectra of high redshift star-forming 
galaxies. 

Subject headings: galaxies: starburst — galaxies: halos — UV: galaxies — 
galaxies: individual (Haroll) galaxies: individual (VV114) galaxies: individual 
(NGC1140) galaxies: individual (SBS0335-052) galaxies: individual (TOL0440- 
381) galaxies: individual (NGC1705) galaxies: individual (NGC1741) galaxies: 
individual (IZW18) galaxies: individual (NGC3310) galaxies: individual (Haro3) 
galaxies: individual (NGC3690) galaxies: individual (NGC4214) galaxies: in- 
dividual (MRK54) galaxies: individual (IRAS19245-4140) galaxies: individual 
(NGC7673) galaxies: individual (NGC7714) 



Introduction 



While ACDM simulations of a universe dominated by dark matter and dark energy 
have led to a robust understanding of the large scale structure of the universe, many 
questions remain when studying smaller scales where baryonic physic s become important 



([Robertson et al.l l2005l : ISommer-Larsen et al.l Il999l : iBenson et al 



20031 ) . Thus, understand- 



ing the interplay between a galaxy's components (e.g. gas, stars, radiation, black holes) and 
its environment is an important next step for cosmology studies. 

In particular, feedback on the surrounding gas-phase baryons provided by massive stars 
is believed to play a major role in the formation and evolution of galaxies. Starburst driven 
galactic winds are the most obvious manifestation of these processes. These "superwinds" 
appear to be omnip resent in galaxi es with star formation rates (SFR) per unit area > 
10~^ Mq yr^^ kpc~^ (jHeckmanll2003l ) . Locally, these starbursts alone account for roughly 20% 
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of th e high-mass star formation and 10% of the radiant energy production flHeckman et al. 
20051 ). Powered by supernovae and stellar winds these outflows drive gas, d ust, energy, and 



meta l s out of galaxies an d possibly into the inter-galactic medium (IGM, iHeckman et al. 



19981: 



Veilleux et al l boosh. The higher overa ll star-formation rates at z ~ 1 and above 



Bunker et aPbooi Iflopkins fc Beacom I [2OO6I ) imply that star-formation driven winds are 
even more important in the early universe. Such winds at high redshift almost certainly 
had a profound imp act on the energetics, ionization, and chemical composition of the IGM 
(lAguirre et al.lboosh . 



Lyman break galaxies (LBGs. lSteidel et al.lll999l ) are the most widely studied population 
of high redshift star forming galaxies as they can be efficiently selected using the Lyman- 
break dropout technique. Rest frame UV spectroscopy of LBGs have shown that starburs t 
driven galactic winds appear to be a ubiquitous property of LBGs ( IShapley et al.l l2003l ) . 
During the cosmological important epoch of 2 ~ 2 — 4 LBGs constit ute a significarit (an d 
possibly dominant) fraction of the population of star forming galaxies (jPeacock et al.ll2000l ). 



Stu dies of the outflows in star forming galaxies are complicated by their multiphase 
nature (jStrickland et al.l |200J) . The mechanical energy from t he supernovae and stella r 
winds create an over-pressurized cavity of hot gas at ~ 10^ K ( IChevalier fc Cleggi Il985l ). 
The hot gas expands and cools adiabatically as it accelerates into lower pressure regions. 
Shock fronts form at the int erface between the hot outflowing gas and the ambient material 
( IStrickland &: Stevens! I2OOOI ). Some of this ambient material is subsequently swept up into 
the outflow itself. Unfo rtunately, directly studying t he hot wind fluid is only possible in a 
few local galaxies (e.g. [Strickland fc HeckmanI 120071 ) due to the low emissivity of the hot 
gas and observational constraints in the hard X-ray band. Observations then of superwinds 
focus on either material entrained in the wind fluid (e.g. neutral, ionized, molecular gas 
and dust) or created t hrough the interface between the hot wind and the ambient material 
(iMarcolini et allboosh . 



The coronal gas (T ~ 10^ — 10^ K) observed in galactic winds is of particular importance. 
This phase is intimately connected to the mechanical/thermal energy that drives the outflows 
and places important limits on the cooling efficiency of the outflow. The coronal phase is 
probed through observations of the O VIAA1032, 1038 doublet found in the rest frame FUV. 
Unfortunately, observations of this wavelength regime in LBGs are generally limited by its 
location deep within the Lyman a forest. Our knowledge of the properties of the winds in 
LBGs is generally then limited to what can be inferred from the interstellar UV absorption 
lines longward of Lyman a. 



Studies of local star-forming galaxies in this wavelength regime have been limited by 
UV absorption in the Earth's atmosphere. Hopkins ultraviolet telescope (HUT) observations 
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of star forming galaxies analyzed by iLeitherer et al. ( 2002) pro bed the 912-1800 A regime 



but with very low spectral resolution. iHeckman et al.l (119981 ) studied a large sample of 



45 starburst galaxies using the International Ultraviolet Explorer (lUE) with similarly low 
resolution and wavelength coverage only longwa rd of about 1150 A. HST UV (> 1215 A) 



spect ral studies have also been undertaken (e.g. IVazquez fc Leithererl l2005l : ISchwartz et al. 



20061 ) although they focus primarily on specific star clusters given the small spectroscopic 
apertures of HST. 

The situation changed however with the launch of the Far-ultraviolet Spectrographic 
Explorer (FUSE) in 1999. FUSE provides high resolution spectroscopy with an aperture 
that is well matched to observing global FUV properties of local star-forming galaxies. Two 
previous w orks using overla pping data sets have used FUSE to study dust at tenuation in 
the FUV JBuat et all 120021 ) and stellar populations JPellerin fc RobertI booj). However 



while s everal previous works have studied specific galaxies with FUSE (e.g. 



2001a: 


Aloisi et al. 


2003 


Grimes et al. 


2006), a St 



2004J : iGrimes et al.ll2005l : iBergvall et al 



Heckman et al. 



2006 



of star-forming galaxies has not previously been undertaken. 

In this paper we will focus on the absorption lines produced by the interstellar medium 
(ISM) and observed in the FUSE wavebands. We will first introduce our sample selection 
method and analysis techniques in Section [2l Next we will discuss the outflow velocity 
distributions of the various ISM absorption features in Section 13.11 Section 13.21 focuses on 
the different physical origin of the coronal gas as traced by O VI. We then compare the 
outflow velocities and strengths of the absorption features to the physical properties of the 
galaxies in Sections 13.31 and 13.41 respectively. Section 13.51 discusses the possible escape of 
Lyman continuum emission in local starburst galaxies. 



2. FUSE Observations 

The FUSE satellite is a high resolution spectrometer with sensitivity in the 905 - 1187 
A range. Four co-aligned prime focus telescopes, each with its own Rowland spectrograph, 
comprise the four FUSE channels. Each channel has its own focal plane assembly with 
four apertures. Only two science apertures are relevant to this work, the low-resolution 
aperture (LWRS; 30" x 30") and the medium-resolution aperture (MDRS: 4" x 20"). Spectral 
resolution (full width at half maximum) ranges from about 17km s~^ for a point source to 
about 100km s~^ for a uniform source filling the LWRS. Two channels have an Al-|-LiF optical 
coating while the other two are SiC providing sensitivity for the ~990 - 1187 A and ~905 
- 1070 A ranges respectively. The spectra are imaged onto two micro-channel plate (MCP) 
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detectors (1 & 2) each with two sides (A & B). The four channels and two sides to each 
MCP produce eight individual spectra for every FUSE observation. Telescope acquisition 
and guiding is controlled by a fine error sensor which is aligned with the LiF 1 channel for 
the observations in this paper. 

2.1. Sample Selection 

One of the primary motivations for this work is to study the properties of O VI absorp- 
tion in star-forming galaxies. The O VI feature is actually a doublet at 1032 and 1038 A. In 
our galaxies the feature is optically thin, and the A1038 line will be half as strong as the A1032 
line. Since the weaker A1038 line is severely blended with C IIA1036 and O IA1039, we do not 
attempt to measure it in our spectra. Therefore, our sample selection is driven by the desire 
to pick star-forming galaxies in which we cow/ii detect O VIA1032. While O VIA1032 is gener- 
ally a prominent feature in these galaxies, it can be very broad and is found in a complicated 
region of the spectra. The wavelength region around O VIA 1032 is particularly intricate due 
to the presence of several other strong features most noteably Ly/5 A1025 and C IIA1036. 
More problematically. Milky Way (MW) O VI and C II are strong absorption features that 
may also fall on the observed redshifted extra-galactic O VIA1032 feature. Therefore, we 
cut galaxies with systemic velocities Vsys < 600kms~^ and 1000 < Vsys < 1500 km s~^ as ex- 
tragalactic O VI would then be blended with the strong Milky Way O VI and C II features 
respectively. The two highest redshift galaxies in the original sample were also dropped as 
the O VI was shifted onto the gap between the LiF channels at ~1090A. We also require 
S/N > 4 in 0.078 A bins in the LiF lA galaxy spectrum. The LiF lA spectrum is the 
highest S/N spectrum of the eight spectra that are obtained during every FUSE exposure. 
This reduced our sample from the 54 star-forming galaxies we found in the archive to only 
15 objects. We have however made an exception for one galaxy, NGC 4214. One of the 
side-effects of dropping galaxies with Vgys < 600kms~^ is that we dropped a majority of 
the lower SFR galaxies found in the archive. In the case of NGC 4214, it is a particularly 
high quality dataset (S/N=14.6) with sufficiently narrow absorption features so that they 
are cleanly separated from the MW features. We therefore include NGC 4214 in our sample 
for a total of 16 galaxies. Table [1] lists the galaxies in our sample along with many of their 
physical properties. 

The FUSE apertures for these observations are generally pointed directly at the location 
of the most intense SF in every galaxy. As the FUV emission of most of the starbursts in the 
sample is strongly peaked near the galaxy center, the LWRS aperture includes the majority 
of the FUV emission for each galaxy. We are therefore probing the FUV properties of the 
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entire galaxy, similar to observations of high redshift LBGs. The lone MDRS pointing, 
NGC 7714 has a compact nuclear starburst so that the majority of FUV emission of the 
galaxy is enclosed even in the smaller aperture. We have overlaid the F?7iS'£' aperture and 
orientation on FUV or UV images for 15 of the galaxies in this sample in Figures [2] and [31 
Tol 0440-381 has been excluded as a UV image of that object is not available. All of the 
images in this sample are continuum UV images although the majority are at significantly 
longer wavelengths than our FUSE spectroscopy. Most of the images are HST WFPC using 
either the F255W (A^ = 2599 A), F300W (A^ = 2987 A), or F336W (A^ = 3344 A) filters. All 
but three of the images are from high spatial resolution instruments with the exception 
of the low resolution images of Mrk54 (SDSS U-Band), NGC 4214 {GALEX NUV), and 
Haro3 [GALEX ^\JV). A logarithmic image scaling has been used to highlight the various 
features of these galaxies although it does somewhat misrepresent the relative strength of 
the central features. Specifically, only NGC 3310 has significant extended emission outside 
of the central regions of the aperture. 

The selected galaxies in our sample represent a wide range in galaxy types and proper- 
ties. The galaxies span over two orders of magnitude in stellar mass and nearly three orders 
of magnitude in SFR. We have included several dwarf galaxies (e.g. IZwl8, NGC 4214, 
NGC 1705), spirals (NGC 3310, NGC 7714), compact ultraviolet luminous galaxies (UVLG; 
VV114, Haroll, Mrk54), and highly IR-luminous mergers (NGC 3690, VV114). The 
UVLGs are of particular interest due to their strong si milarities to LBGs in terms of their 



masses, UV luminosities, and UV surface brightnesses (IHeckman et al.ll2005l : iHoopes et al. 



20061 : lOverzier et al.ll2008[ ) 



2.2. Data Reduction 

We have chosen to use a consistent data reduction method throughout the sample aimed 
at maximizing S/N. This method involves several extra analysis steps but provides superior 
data quality compared to the standard data reduction method for the targets in our sample. 
While some of the observations have high S/N ratios in the LiF channels, the lower effective 
area SiC channels are almost always relatively low S/N. Our study focuses principally on the 
LiF channel data but we use the SiC data for measuring the very strong C IIIA977 absorption 
line and investigating possible escaping Lyman continuum emission. 

The first step was to run all raw exposure data through the CalFUSE 3.1.8 data pipeline 
to produce intermediate data (I DF) files. We t hen ex tracted all eight flux calibrated spectra 



for every FUSE exposure. See iDixon et al.l (120071 ) for a complete discussion of the cali- 



brations and corrections applied to the data. Mirror misalignment and thermally-induced 
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rotations of the spectrograph gratings cause small zero point shifts in the wavelength calibra- 
tion. Therefore we cross calibrate the flux calibrated exposures to determine the individual 
zero-point wavelength shifts. Given that most of our observations are relatively low S/N 
spectra of extended sources it is difficult to accurately determine many of these wavelength 
offsets. We subsequently have chosen conservative offsets. A histogram of the offsets used 
for all of the individual exposures in shown in Figure [H The majority of the offsets are 
< 0.013 A (< 4kms-i) which is significantly less than the wavelength inaccuracies caused 
by localized distortions of the detector (0.025 A ~ 8kms~^). 

After determining the wavelength offset for each individual spectra we retroactively 
apply the correction to the individual IDF files. We also examine the count rates of the 
individual exposures and throw out the exposures with significant discrepancies. These 
variations are usually caused by misalignment of a mirror on the target. The IDF files for each 
exposure are then combined to create one total IDF file for each channel. We then re-extract 
a combined fiux calibrated spectrum which also has the zero-point wavelength corrections 
applied. More importantly, however, the longer, combined IDF files allow the CalFUSE 
pipeline to do an improved background subtraction. In particular, the pipeline fits the 
various backgrounds (e.g. detector, airglow, and geo-coronal scattered light). For the shorter, 
individual exposures it simply scales the background files by exposure time when extracting 
fiux calibrated spectra. As there can be significant background variations, exposure time 
scaling can introduce significant background subtraction errors. The backgrounds dominate 
the signal for many of these spectra so that errors in the background subtraction introduces 
systematic offsets in the flux scale. Even when using background fltting, these systematic 
offsets can still be seen in many of the low S/N SiC spectra. 

In order to minimize the detector background, we ha ve extracted our spectra using the 



smaller, point source detector regions (IDixon et al.l 120071 ) for all objects in our sample with 



the exception of NGC1705, NGC3310, and NGC4214. A visual examination of the 2-D 
detector images shows that these three galaxies have source counts falling outside of the 
point source aperture. These are the three highest flux galaxies within our sample. Figure 
[2] suggests that at least for NGC3310, the galaxy has considerable extent within the FUSE 
aperture. By using the smaller point source region for the remaining 14 galaxies we include 
the same amount of source flux but minimize the included background. 

For galaxies which contain a signiflcant fraction of night only data, we have thrown out 
the daytime data. Signiflcant airglow features contaminate the day time spectra and are 
best avoided whenever possible. Table [2] lists the FUSE observational data for all of the 
galaxies included in our sample. 
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2.3. Data 

The SiC 2A, LiF lA, and LiF 2A spectra for every galaxy in our sample are shown in 
Figures H] - [191 These channels cover most of the wavelength regime between 920 - 1180 A. 
We use LiF 2B for these panels instead of LiF IB as it is rarely affected by the " worm". The 



"wor m" is depression in the observed flux that can span over 50 A in length ( Dixon et al. 



20071 ). Shadowing by a grid of quantum efficiency wires where the secondary focus falls above 
the detectors and close to the location of the wires causes these observed drops in flux. This 
feature is seen in almost all of the LiF IB spectra in our sample but does not affect any of 
the other channels. 

Most of the common MW, galactic, and airglow features are marked on these figures. 
The galactic absorption lines probe the stellar population, interstellar medium (ISM), and 
outflowing gas. H2, usually of MW origin, is present in many of the galaxies and is a 
significant hurdle in identifying and fitting many of the weaker ISM absorption features. 



2.3.1. Gaussian Fits 
We have fit absorptio n lines to mo st of the rn ost prominent features using the IRAF 



(ITody fc Fitzpatrickl Il996l ) tool specfit (jKrisd [l99J) . We have not attempted to fit every 



possible feature but have instead focused on the strongest, separable lines that we positively 
identify (this can be particularly difficult in galaxies contaminated by MW H2). Each line was 
fit using a freely varying powerlaw for the continuum and a symmetric gaussian absorption 
line. When absorption lines are blended (e.g. O I A989 and N III A990) we fit them 
simultaneously. Due to the differences in the redshifts of the galaxies (and hence differing 
degrees of blending or contamination by lines due to the Milky Way), it was not possible to 
consistently define the spectral windows over which the adjacent continuum was fit for any 
given absorption feature. 

We note that the compact UV sizes of most of these galaxies implies that the velocity 
resolution in our spectra willi typically be roughly 20km s~^ (FWHM). The worst resolution 
will be for NGC 3310, which based on its angular size we estimate will be about 60km s^^. 
In all cases, the observed absorption-lines are well-resolved. 

For the majority of absorption lines we have used the LiF channels to make two in- 
dependent measurements of the equivalent width, velocity, and full-width half-maximum 
(FWHM). For lines with only one measurement, the line either fell in the gap between LiF 
channels or was in the region (1125-1160A) affected by the worm on the LiF IB channel. We 
have generally ignored the SiC channels above 1000 A due to their lower S/N. The results 
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of our fits are listed in Tables [3] - [151 



The listed statistical errors are la and are calculated from the minimization error matrix 
which is re-scaled by the reduced value (iKrisd Il994l ). Systematic errors are probably 



not large in the FWHM measurement but may be important to the equivalent width and 
centroid velocity. The equivalent width determination is largely a problem for saturated 
lines in low S/N spectra. Background over-subtraction is a common problem for these 
spectra so it is difficult to determine the continuum level and fiux zero point. Systematic 
wavelength calibration errors exist and are pa r ticula rly hard to quantify for our extended 
sources using the LWRS aperture. iDixon et al.l (120071 ) suggest that they are order 10 kms~^ 
which seems consistent with the observed variations in our results. A comparison of the 
derived parameters for features with data on two LiF channels shows broad agreement, 
generally well within the la error bars as expected. In these cases, we have used the error 
weighted value to create a single measurement for the analysis that follows. 

Close up plots of six strong absorption features for all 16 galaxies are shown in Figures 
I2U] - [33 The plots focus on the ISM absorption lines observed in almost all of the spectra, 
specifically O VIA1032, C IIA1036, C IIIA977, O IA989, N IIA1084, and O IA1039. We have 
displayed these spectra even when the absorption feature was not observed. MW ISM 
contamination, H2, airglow, instrument coverage gaps, and signal quality play an important 
role in which features are actually detected. Line blending is also important most particularly 
for O IA989 as it falls right next to the strong N IIIA990 feature. 

Another potentially problematic blending is C IIIA977 with O IA976. However, we 
have checked the significance of the contamination by modeling the O IA976 line using 
our gaussian fits to O IA1039. From iMortonI (Il99ll ) we expect O IA976 to be a factor of 
2.7 weaker than O IA1039 for an optically thin feature. In Figure [Ml we have overlayed an 
O IA9 76 absorption feature with the same velocity and FWHM as the O IA1039 profile but 
we have lowered the equivalent width by a factor of 2.7 for NGC 1140 and NGC 3310. As 
shown in these plots, this feature is significantly weaker than the observed C IIIA977 profile. 
Fitting these profile with the O IA976 contribution does not change the derived equivalent 
widths, centroid velocities, and FWHMs in a significant manner. The largest percentage 
change is seen in the FWHMs which lower by 12 and 11 kms~^ for NGC 1140 and NGC 
3310 respectively. This is only 25% of the la statistical error. We therefore conclude that 
the weak O IA976 is not significantly affecting our fits to the C IIIA977 profile. 

The stren gth of the absorption lines varies widely across the sample. As discussed in our 
previous work [Grimes et al.l (|2006[ ) and [Grimes et al.[ (120071 ) , most of the strong absorption 
line profiles are non-gaussian. The absorption profiles frequently show significant structure 
suggesting that in many cases these are an aggregate profile representing multiple absorbers 
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spread throughout the galaxy's ISM and starburst driven outflow. Many of the profiles are 
asymmetric, with stronger absorption on the blueward side of the line core. This is consistent 
with absorbing gas entrained in and accelerated by an outflowing wind fluid. 

In our previous work we have successfully modeled many of the stronger lines with two 
gaussian component fits. Physically, for both Haro 11 and VV 114 the stronger gaussian was 
associated with the host galaxy ISM while the blueshifted gaussian represented outflowing 
gas at several hundred kms~^. For this work, covering a wide range of galaxy star formation 
rates and data quality, a multiple component fit is not feasible. Therefore, we have used two 
additional methods to study the strongest absorbers in our sample. 



2.3.2. Apparent Optical Depth Method 



For O VIA1032 we use the apparent optical depth method as described in lSembach et al. 



( 120031 ). This method assumes that the O VI absorption features are not saturated - which 
is consistent with the observed O VI profiles. The apparent column density as a function of 
velocity is given by 

N.(.) = = 3.768 X 10-^ (1) 

^ ^ /A 7re2 /A(A) ^ ^ 

where 

ra(z/) = -ln[IobsM/IoM]. (2) 

lo(z^) is the continuum intensity while lohs{^) is the observed line intensity. The O VI line 
profiles in Figures [20] - [351 show that the line is clearly resolved when present. Therefore we 
can obtain the column density by integrating over velocities as N = jN(z/)dz/ = jNa(z/)dz/. 

The apparent optical depth method also allows us to calculate the velocity centroid and line 

1 

width using z/ = / vTs^{i')dv / J rg^{v)di' and 6 = [2 /(z/ — p)^ra(z/)(iz// / Tg^{u)du\'2 respectively. 
Table [19] lists our results. 



In our previous work we have calculated the O VI parameters for VV114 and Haro 11 
from gaussian fits. For VV 114 these two methods are in excellent agreement as we had found 
No VI ~ 2.0 X lO^^cm"^ while we now calculate Nqvi ~ 2.2 x lO^^cm"^ using the apparent 
optical depth measurement. These results are consistent because the O VI profile shape is 
approximately gaussian for VV 114 (see Figure [2T]) . We find however a larger discrepancy 
for Haro 11 with Nqvi ~ 4.2 x lO^^cm"^ and Nqvi ~ 8.0 x lO^^cm"^ for the gaussian and 
apparent optical depth methods respectively. For Haro 11 (Figure [20]) . a multicomponent 
profile is observed. In this case the apparent optical depth method better includes the ex- 
tended blue wing. We also now derive a larger absorption line width and a more blueshifted 
velocity centroid, consistent with a significant contribution from the extended blue wing of 
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the absorption profile. The apparent optical depth method provides a significant improve- 
ment over gaussian profiling for O VI and is applicable as O VI is not saturated or blended 
significantly with other lines (partly due to our sample selection). Table [T^ does not include 
four galaxies. For NGC 1140, 1 Zw 18, and IRAS 19245-4140 we do not positively detect O VI 
absorption. In Tol 0440-381 the O VI is a broad, shallow feature making it impossible to 
reliably estimate the velocity limits over which the above integral should be performed. 



2.3.3. Non-parametric Equivalent Widths 

Several other strong absorbers exist in the spectra of these star-forming galaxies, specif- 
ically C IIA1036, C IIIA977, and N IIA1084. Unlike O VIA1032 however, these lines are 
frequently saturated at line center and therefore we are unable to apply the apparent optical 
depth method. We th erefore measure the equiva lent widths and velocity centroids of these 



lines using the IRAF (ITody fc Fitzpatricklll996l ) tool splot. The procedure subtracts a user 



defined linear continuum and sums the resulting flux to calculate the equivalent width. The 
comparably strong O IA989 and N IIIA990 lines are frequently blended, so these lines have 
been excluded from the procedure. Our results are listed in Table [201 These lines are fre- 
quently saturated, particularly C III and C II, so the equivalent widths are essentially just 
a measure of how broad the absorption feature is. 

We have also compared the non-parametric equivalent widths, centroids, and widths 
with those derived from our previous gaussian fits in Section [2.3. 1[ There is broad agreement 
for the derived values using the C II and N II features. The principal outliers are the C II 
centroid velocities for some of the galaxies with higher SFR (> a few 10^° Lq) where the 
outflow velocities are found to be higher from the non-parametric fits. In the case of C III, 
there is a systematic shift to higher outflow velocities for the non-parametric values. This 
is observed in almost all of the galaxies in our sample not just in the galaxies with higher 
SFR. These findings are consistent with our motivation for using a non-parametric method 
to measure these absorption features. Specifically, as line strengths increase the absorption 
features are found to be increasingly non-gaussian and blue-asymmetric. 

For the analysis that follows unless otherwise stated, we use these non-parametric mea- 
surements for the C III, C II, and N II features. For O VI we use the results derived from 
the apparent optical depth method and for the remaining lines we use the gaussian data. 
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3. Discussion 

This paper is intended to be a general study of the FUV properties of starburst galaxies. 
As stated previously, we will focus on the strongest ISM absorption features in order to 
include the majority of galaxies in our sample. In particular, we will focus on the following 
lines: O VIA1032, C IIIA977, S IIIA1012, N IIA1084, C IIA1036, O IA1039. These lines are 
detected and measured throughout most of our sample. They also represent three different 
phases of the ISM - coronal (O VI), HII (C III, S III, N II) and HI (C II, O I). 

3.1. Outflow Velocity Distributions 

In Figure EZI we have plotted histograms of the relative velocities for the strongest of 
these absorption features. These lines are listed in order of increasing ionization potential. 
With the exception of the neutral oxygen ion, the other absorption features show the presence 
of outflows ranging up to a few hundred kms~^. As the ionization state increases the lines 
generally have larger outflow velocities. The sole exception is C II which is mildly more 
blueshifted than N II. The C III and O VI lines however are significantly more blueshifted 
than the two lower ionization state ions. 

In Figure [38] we have directly compared the velocities of these absorption features. 
The left panel shows C II versus N II outflow velocity while the right panel is C III versus 
O VI. For the majority of galaxies in our sample, C II is clearly more blueshifted than N II. 
Similarly, O VI has an higher outflow velocity than C III. This confirms the pattern seen in 
Figure [371 

As we have previously discussed, the O VI line (unlike C III, C II, and N II) is not 
saturated for any of the galaxies in our sample. Moreover, an examination of the O VI profile 
shape suggests a different origin for this feature. The strongest part of the O VI profile is 
blueshifted, and the absorption at or near the galaxy systemic velocity is weaker. In contrast, 
the C III, C II, and N II lines are broad with extended blueshifted wings, and have the 
strongest absorption (the line core) at or near the galaxy's systemic velocity. Furthermore, 
the maxiumum outflow speed (defined by the blueward extent of the absorption-line profile) 
is similar for O VI and C III. 

Thus, the gas traced by the O VI line appears to be primarily in the outflow rather 
than in the quiescent ISM, while the lower ionization lines are produced by both. This 
difference is not surprising since O VI has a much higher ionization potential and is produced 
in significantly higher temperature gas. We will discuss this further in Section 13.21 
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Several previous studies have suggested that the me asured outflow v elocit ies increase 
with the ionization potential of the species in starbursts. IVazquez et al.l (120041 ) used HST 
STIS spectroscopy of NGC 1705 (also represented in our FUSE sarnple) a nd found that 
several high ionization lines had higher outflow speeds. ISchwartz et al.l (120061 ) similarly used 
HST STIS to observe starclusters in 17 star-forming galaxies (including NGC 3310 and NGC 
1741 also in this work). While C IIA1335, Al IIA1671, and Si IIA1260 all trace the neutral 
hydrogen phase, they found evidence that C IIA1335 had higher outflow speeds than the 
two slightly lower ionization features. We als o noticed similar b ehavior previously when 
studying the FUSE observatio ns of Haro 1 1 (j Grimes et al.l |2006| . Haro 11 is also in this 



sample). ISchwartz et al.l (120061 ) suggest that this could be explained by less dense, hotter 
gas escaping the starburst region and expanding at a faster speed than the denser, lower 
ionization gas. 

This picture a ppears to be d i fferen t than what is observed in high redshift Lyman break 
galaxies (LBGs). iPettini et al.l (120021 ) have examined high quality spectra of the lensed 
LBG, MS1512-cB58. They found the same velocity distri butions for bot h the l ow and high 
ionization features. Using a composite LBG spectrum, IShapley et al.l (120031 ) also found 
similar outflow velocities for both the lower and higher ionized absorption features. 

As discussed previously C III , with an ionization potential of 24.4 eV, has signiflcantly 
higher outflow velocities than C II (11.3 eV), N II (14.5 eV), and O I. This is then consistent 
with the idea that higher ionization lines have higher outflow velocities. On the other hand, 
the lower ionization C II has a slightly higher outflow velocity than N II (e.g. left panel of 
Figure [38]) . To study this further we have plotted the C IIA1036 versus S IIIA1012 velocities 
(Figure 139!) . S III has a much higher ionization potential (23.3 eV) than C II but is a weaker 
feature. It is not saturated and is well flt by a gaussian model (unlike the stronger features). 
C II, C III and N II all have signiflcantly higher outflow velocities than S III , which is not 
consistent with an increase of velocity with ionization energy. 

In Tables |3] - [18] we listed the absorption features in order o f decr easing values of 
log{XfN/NH)- The oscillator strength (/) values are from iMortonI (|l99ll ) and solar rela- 
tive abundances are assumed for N/Nh- While the majority of the galaxies in our sample do 
not have solar metallicity (Tabled]), this is still a useful predictor of the relative strengths 
of the absorption features. With the exception of O VI (which we have already argued has a 
different physical origin than the other ISM absorption features) the values of the equivalent 
widths follow decreasing log^XfN/Nfj) (see also Table [20]). It is particularly striking that our 
previous ordering of outflow velocities follows the same general trend as the line strengths. 
Both outflow velocities and line strengths decrease from C III, C II, N II, to O IA1039. 



An examination of the line proflles show that the stronger lines have extended blue 
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absorption wings. In fact, these extended absorption wings were our prime motivation in 
using a non-parametric method for measuring the equivalents widths and centroid velocities 
(Section 12.3.31) of the strongest features. This was not necessary for the weaker lines like 
O IA1039 and S IIIA1012 which are well fit by gaussian modeling. As the absorption lines 
increase in strength, the extended blue wing in the absorption-line profile becomes more and 
more conspicuous. The simplest explanation for this is that the high velocity outflowing 
gas has a lower column density than the more quiescent gas in the ISM. We therefore only 
detect the fastest moving gas in the strongest absorption features. This explanation is 
consistent with the observa tions of the LBG cB 58 where the stronger lines span progressively 
broader ranges in velocity (jPettini et al.ll2002l ). It also explains the high er outflow speeds for 
C IIA1335 relative to Al IIA1671 and Si I IA1260 JSchwartz et all booeh as C IIA1335 is the 
strongest feature. Lastly, in the work by IVazquez et al.l (120041 ). the lines with the highest 
measured equivalent widths are preferentially those with the highest ionization energies. This 
result thus appears to be consistent throughout a wide range of galaxy types and absorption 
features. 

This result is particularly relevant to outflow velocities derived in other wavelength 
regimes. For example, the Na IAA5890, 5896 doublet has log{XfN/NH) = —1.94, similar to 
N IIA1084. This suggests that observations of local outflows using Na I absorption are not 
biased by measuring lower ionic states than the FUV and UV measurements of higher ioniza- 
tion absorption lines. The derived outflow velocities should be consistent with those of other 
similarly strong absorption features in the neutral and ionized gas. This could strengthen 
studies comparing observations of Na I absorption in local galaxies with the properties of 
outflows observed in LBGs. Unfortunately, there is only on e galaxy in our sarnple ha ving 
published data on the Na I absorption feature: NGC 4214 (jSchwartz fc Martini |2004| ). In 
this case the agreement is good: they flnd an outflow speed of 23 kms~^, while we flnd 
speeds ranging from to 58 kms~^ depending on the particular line. 



3.2. Origin of O VI Absorption 

The O VI line is a particularly interesting one due to its potential importance in radiative 
cooling and the relatively narrow temperature range over which it is a signiflcant (e.g. T ~ 
10^ — 10^ K). As we have discussed, the kinematics of the coronal gas as traced by the 
O VI absorption line are different from that seen in the cooler gas. In particular, O VI is 
the most blueshifted ion and (on kinematic grounds) is located primarily in the outflowing 
gas and not in galaxy itself. 



Several previous papers have studied O VI in speciflc galaxies, e.g. NGC 1705 (iHeckman et al 
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2001ah . VV 114 flGrimes et alJbood ). and Haro 11 flGrimes et alJboOTf ). They attributed the 



production of O VI to the intermediate temperature regions created by the hydrodynamical 
interaction between hot outrushing gas and the cooler denser material seen in Ha images. 
Such a situation is predicted to be created as an overpressured superbubble accelerates and 
then fragments as it expands out of the galaxy and again as the wind encounters ambient 
clouds in the galaxy halo. 



Heckman et al.l (120021 ) derived a simple and general relationship between the O VI col- 



umn density and absorption line width which will hold whenever there is a radiatively cool- 
ing gas flow passing through the coronal temperature regime. They showed that this simple 
model accounted for the properties of O VI absorption line systems as diverse as clouds in 
the disk and halo of the Milky Way, high velocity clouds, the Magellanic Clouds, starburst 
outflows, and some of the clouds in the IGM. Using the the va l ues fo und in Table dni in 
Figure Ho] we have updated the original plot from lHeckman et al.l ( 12002| ). Our results for the 
O VI measured for the galaxies in our sample are fully consistent with the model predic- 
tions. The fact that starbursts tend to lie at slightly lower column density than the model 
predicts is not surprising, given that the model is so generic and idealized (the exact value of 
the coefficient relating column density to flow speed depends weakly on the exact nature of 
the flow, as described in Heckman et al. (2002). We conclude that O VI traces radiatively 
cooling gas produced in an interaction between the hot outrushing gas seen in X-rays and 
the cool clouds seen in Ha emission, and that this is is a generic feature of starburst galaxies. 



In the case of NGC 1705 iHeckman et al.l (l2001al ) showed that the overall mass of the 
coronal phase gas was signiflcantly smaller than the lower ionization material (by about two 
orders-of-magnitude). They showed that this is c onsistent with the pictu re in which this 
gas is rapidly cooling (it is transitory). Following iHeckman et al.l (l2001al ) we can make a 
rough estimate of the gas mass for our sample of galaxies. This mass is simply the product 
of the O VI column density, the surface area of the O VI region, and the inverse ratio of 
O VI to total column. The most physically plausible idealized geometrical model is one 
in which the coronal gas traces the interface b etwee n spherically expanding wind and the 
ambient gas in the galaxy halo. iGrimes et al.l (120051 ) found that the wind radius as traced 
in soft X-ray emission increased systematically with both the star formation rate and galaxy 
stellar mass. Over the range spanned by our sample, the radius would range from about 0.7 
kpc to about 10 kpc. The conversion from O VI to total column density would be factor 
of about 10^ for gas in collisional ionization equilibrium at T = 3 x 10^ K and with a solar 
oxygen abundance. The implied masses of coronal gas range from about IO^Mq for the dwarf 
starbursts like NGC 1705, NGC 4214, and Haro 3 to about IO^Mq for the most powerful 
starbursts like VV 114 and NGC 3690. 
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It is also interesting to note the lack of O VI emission in the sanaple. While weak 
O VI emission is observed in Haroll (iBergvall et al.ll2006l : iGrimes et al.l 120071 ) there is no 
convincing evidence for it in any of the other galaxies in this sample. Very broad (few 
thousand km s~^) O VI emission features appear in a few of t he spectra, consistent with 
stellar winds in the population of O stars (jPellerin et al.ll2002l ). As O VI emission is the 
dominant cooling mechanism for coronal gas, the general lack of observed emission implies 
that radiative cooling is not significantly i mpedi ng the growth of the outflow. This is true 
even in the case of Haroll iGrimes et al.l (120071 ). They suggested that the O VI emission 
may not be from coronal gas cooling but instead could be resonant scattering of far-UV 
continuum photons from the starburst off O VI ions in the backside of the galactic wind 
traced by the blueshifted O VI absorption. 



3.3. Dependence of Outflow Velocities on SFR &; Stellar Mass 



In IGrimes et al.l ( l2005l ) we analyzed Chandra data of star-forming galaxies ranging 
from dwarf starbursts to ultra- luminous infrared galaxies (ULIRGs). We showed that the 
temperature of the hot gas was not dependent upon the galaxy's SFR. The thermalized X-ray 
gas is thought to be pr oduced in a reverse shock b etween the hot out-rushing wind fluid and 
the ambient material (IStrickland fc StevensI l2000l ) . and its temperature is then determined 
by the wind outflow velocity. Since a relationship between X-ray temperature and SFR was 
not observed, it follows then that the wind outflow velocity does not depend on the SFR. 



In contrast, iMartinI (120051 ). ISchwartz et al.l (120061 ) . and lRupke et al.l (120051 ) hav e shown 
that the o utflow speed r neasu red in absorption lines does depend on the SFR. iMartin 
(120051 ) and iRupke et al.l (120051 ) observed outflows in luminous infrared galaxies (LIRGs) 
and ULIRGs using the Nal D doublet and found that the outflow velocities increased with 
the SFR up to values of SFR of a few tens yr~^ At higher SFRs (extending up to sev- 
eral hundred Mq yr~^), the relation appears to saturate with the outflow velo cities no longer 



increasing. Similarly, in the STIS starburst sample discussed in Section l3.1llSchwartz et al. 
(120061 ) also found that the outflow speeds increased with the SFR below ~10 M0 yr~^. 



Martini (120051 ) suggests that the simplest explanation is that there is not enough mo- 



mentum carried by the wind fluid in the lower SFR galaxies, so that the winds are unable to 
accelerate the entrained absorbing material up to the velocity of the wind fluid itself. The 
relation bet ween SFR and the outflow speed seen in the absorption-lines saturates at about 
500 kms~^ ( Martini 2005 ) and this is reached at starburst luminosities sim few xlO^^L0. 
If the X-ray emission in galactic winds arises in wind shocks, the observed te mperatures of 
kTx ~ 0.3 to 0.7 keV correspond to wind shock speeds of ~ 500 to 800 kms~^ (IGrimes et al. 
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20051 ). This is then compatible with the idea that at sufficiently high SFR (high enough 



wind momentum) the entrained absorbing gas has been accelerated to near the wind speed. 

We can check these results with our data. To that end, we will use the sum of the IRAS 
IR and FUSE UV luminosity as our tracer of the SFR (see Table [1]). The UV luminosity 
is based on the continuum flux at 1150 A and is corrected for Milky Way extinction. The 
combined UV and IR luminosity provides a reasonable estimate of the intrinsic UV flux of 
the galaxy, and thus traces recent SF. 



Martini (120051 ) also found that the outflow velocity increased with the rotational speed 
of the galaxy. This suggests that faster outflows are found in more massive galaxies. We 
check this correlation using the near-infrared luminosity of the galaxy (Lk) as a proxy for 
stellar mass. The K-Band Cousin s-Glass- Johnson (CGJ) luminosities are derived from the 



2MASS atlas flJarrett et al.l 120031 ). The 2MASS K and J Ba nd fluxes given in the catalog 



were converted to K-Band CGJ using the transformations from lCarpenterl (120011 ). We expect 
the K-band luminosity to trace the stellar mass to a factor of roughly two (e.g. Leitherer & 
Heckman 1995; Bell & de Jong 2001; Dasyra et al. 2006). 

We have made no attempt to correct for the viewing angle dependence of the outflow 
velocity. This is primarily a practical decision as we do not have a robust way of obtaining 
the orientation of the outflow in these galaxies. While this will increase the scatter in our 
outflow velocity plots, sample selection effects probably minimize the resulting damage. In 
particular, an edge-on spiral galaxy would not have had enough FUV flux to meet our S /N 
requirements as the disk heavily absorbs the nuclear starburst (e.g. M82). Spiral arm 
features are apparent for several of the galaxies in our sample (Figures [2] and [3]) suggesting 
that we are observing the galaxies roughly face-on, and hence roughly along the expected 
outflow direction. Several other of the galaxies have irregular morphologies. In these irregular 
galaxies it is not at all apparent how inclination angle would be defined. As an example, 
Haro 1 1 has several regions of intense SF which high resolution HST ACS imaging suggests 
could be driving multiple outflows. 

In Figures SI] and 112] we have plotted the O VI outflow velocity as a function of the near- 
IR K-band luminosity and the UV+IR luminosity, using two measurements of the outflow 
velocity. The ffist plot uses the median velocity defined by the apparent optical depth 
method (e.g. half the absorbing column density lies on either side of this velocity). The 
second plot traces the blueward extent of the outflow (the 'maximum' outflow speed). We 
define this as that velocity for which 90% (10%) of the absorbing column density lies to 
the red (blue). The maximum outflow velocity correlates better both the K-band and the 
UV+IR luminosity that does the median velocity. For comparison, in Figure US] we have 
plotted the C III outflow velocities (medians) as a function of 2-Mass K-band and UV+IR 
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luminosity. 

The outflow speed of O VI clearly correlates better with with the SFR than does C III. 
The increase in outflow velocity with increasing SFR is also observed in the other strong 
absorption features such as C II and N II. We have noted previously that the profile shapes 
of these strong lines suggest that there are two components (higher column density gas at or 
near the systemic velocity and lower column density outflowing gas). In contrast, the O VI 
absorption seems to originate only in the outflow. This may explain the stronger correlation 
ibetween SFR and outflow speed seen in O VI. 

Since the outflow speeds increase with both SFR and stellar mass, it is instructive to 
examine the outflow speeds as a function of the specific SFR (SFR/M*). To do this, we 
have plotted the median O VI and C III outflow velocities versus (Ljr + Luv)/Lk in Figure 
SH While these plots show a significant amount of scatter, we see that the outflow speeds 
generally increase with the specific SFR. The O VI outflow velocity correlates more strongly 
with (Lin + Luv)/Lk than does the C III velocity. As previously discussed, the properties of 
O VI are closely tied to the wind and thus presumably to the SFR. 

The case of NGC3310 requires further discussion as it is a significant outlier in most of 
our plots, with an outflow speed that is larger than expected for its SFR or mass. We note 
that the very large outflow speeds we see have been confirmed with STIS far-UV spectra of 
three individual regions in NGC 3310 (Schwartz et al. 2006). Examining the UV morphology 
of NGC 3310 in Figure [21 we see that NGC 3310 is the only galaxy with significant star 
format ion spread throughout multiple locations in the FUSE aperture. Radio HI velocity 



maps (IKregel fc Sancisill200ll ) of the regions observed by FUSE display a spread of almost 
200 kms~^ which could increase the observed width of our absorption features, but should 
not produce highly blueshifted absorption. Another way in which NGC 3310 stands out 
within our sample is that it is the highest metallicity galaxy with log(0/H) + 12 = 9.0. In 
principle, higher metallicity could allow us to probe gas with lower total column density. 
However, as we will show below the strengths of the strong absorption lines do not depend 
significantly on metallicty in our sample. Thus, this is unlikely to explain the large observed 
outflow speed in NGC 3310. 

NGC 3310 is inclu ded in the young st e llar p opulation study of 16 starburst galaxies (11 



also in our sample) by iPellerin fc RobertI (120071 ). By fitting the FUV spectra with models 



of instantaneous bursts, they estimated an age of 18 ± 2 Myr for the starburst in NGC 3310 
which was 10 Myrs older than any of the other galaxies in their sample. As the Type II 
supernovae that drive the galactic wind persist for 40 Myr, the total amount of kinetic energy 
injected into the wind will increase as the starburst ages. NGC 3310 appears to be a more 
mature starburst and may therefore have a more evolved wind. 
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3.4. Absorption-Line Equivalent Widths 

It has known for some time that the stren gths of the UV interst eUar absorption-lines 



correlate with other key properties of starbursts. iHeckman et al.l (Il998l ) used lUE spectra of 
45 local starbursts to show that the equivalent widths of the strongest lines arising in the HI 
phase of the ISM correlated postively with the amount of UV extinction and reddening, the 
galaxy rotation speed and absolute magnitude, the metallicity, and SFR. Unfortunately, the 
poor spectral resolution and modest S/N of these spectra did not allow them to determine the 
astrophysics underlying these corr elations. Because sta rbursts in more massive galaxies are 



known to have higher metallicity (ITremonti et al.ll2004l ). higher SFR, and higher extinction 



( iHeckman et al.l Il998l ) it is especially difficult to sort out the fundamental relations from 



secondary correlations. 

In general, the equivalent width of an interstellar absorption-line will depend on very 
different parameters depending on the optical depth of the line. In the limiting case of 
optically-thick gas the equivalent width is set purely by the product of velocity width of 
the line (the Doppler b-parameter) and the fraction of the UV continuum source covered 
by the absorber (the covering factor). In the case of optically thin gas the equivalent width 
is determined by the product of the ionic column density, the oscillator strength of the 
transition, and the covering factor. In our data we have examples of both optically thick 
and optically-thin transitions. Because the optically-thick lines are black (or nearly so) at 
line center, we know the covering factor is near unity. 

We also probe gas spanning a range in ionization potential. We have described the 
properties of the coronal gas traced by O VI at some length, so here we will confine our 
attention to the lines tracing the HI phase (ionization potentials less than a Rydberg) and 
the HII phase (ionization potentials between 1 and 4 Rydbergs). For each of these two 
phase we will use one optically-thick and one optically-thin line. Hence, in Figure HS] we plot 
the equivalent widths of O IA1039, C IIA1036, S IIIA1012, and C IIIA977 versus gas phase 
metallicity, K-band luminosity, and IR+UV luminosity. The first two lines trace the HI 
phase, with the O I line being optically-thin and C II optically-thick. The second two lines 
trace the HII phase, with the S III line being optically-thin and C III optically-thick. 

In fact, the only pattern that emerges is that the poorest correlations are those between 
the equivalent widths of the lines tracing the HII phase (the C III and sthw lines) and 
the metallicity (Kendall = 0.38 and 0.44 respectively). All the other correlations are 
of comparable strength (Kendall Ta = 0.56 to 0.75). This may refiect the relatively small 
sample size and the fact that the different starburst propertiues are interconnected: galaxies 
with higher masses will have higher velocity dispersions, higher SFRs, and higher metallicity. 



- 20 - 



3.5. Lyman Continuum Emission 



One of the major puzzles in cosmology today is determining the source of t he UV 
photons required to reionize the universe at redshifts 6 < z < 14 (IFan et al.l l2006l ). The 
known population of AGN do not ap pear to be sufficient by a wide margin, an d the leading 
candidates are star forming galaxies (IStiavelli et al.ll2004j : iPanagia et al.ll2005l ). While stars 
in the early universe are thought be able to produce enough UV photons, it is not well 
understood what fraction of these photon s can escape through the neutral hydrogen gas of 



the host galaxy (lYan fc Windhorst 



20041 ) into the IGM. Superwinds, driven by starbursts, 



may be able to e nhance the fracti on of ionizing radiation by clearing channels through the 
surrounding gas (IDove et al.ll2000l ). 



The escape o f ioniz ing continuum radiation from star forming galaxies wa s first observed 
bv ISteidel et al.l (120011 ) in a sample of z ~ 3.4 LBGs. More recent work by IShapley et al. 



( 120061 ) focused on deep rest-frame UV spectroscopy and found a sample-averaged relative 
Lyman continuum escape fraction [fesc) of 14%. They however detected Lyman continuum 
emission in only 2 of the 14 galaxies in their sample, which suggests that even at high 
redshifts (z ~ 3) ionizing radiation is either relatively rare or viewing angle dependent. 

Determining the escape fraction in starburst galaxies is therefore an important step in 
evaluating the importance of starburst galaxies in reionizing t he universe. Lo c al sea rches 
for Lyman continuu m emission have so far been inconclusive. iLeitherer et al. I (Il995h and 
Hurwitz et al.l (119971 ) found fesc ^ 10% in a galaxy sample based on i^t/Tobservationsi and 
Deharveng et al. (2 001) obtained a similar result for Mrk 54 using FUSE. Most recently 
Bergvall et al.l (120061 ) published an analysis showing the det ection of Lyni a n coii tinuum emis- 
sion in Haroll. We however re-analyzed the same data in iGrimes et al.l (120071 ) and did not 
find any convincing evidence of Lyman leakage. We showed that uncertainties in the back- 
ground subtraction and geocoronal contamination were larger than any detected signal. We 
therefore placed a limit on the escaping Lyman continuum emission in Haro 11 of fesc ^ 2%. 



3.5.1. Below the Lyman Limit in Other Starbursts 

We have looked directly for Lyman continuum emission in other galaxies observed by 
FUSE. In order to move the Lyman limit farther from the detector edge, we have chosen the 
nine highest redshift galaxies in our sample. Figures |16] and HT] show the SiC 2A spectra for 
the nine galaxies plotted in order of increasing redshift. While the SiC IB spectra also cover 
this wavelength regime, background subtraction problems are even more severe for the SiC 
IB channel. In order to minimize contamination by dayglow we have plotted the night only 
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spectra for all of the galaxies with the exception of VV 114. In the case of VV 114, there was 
not enough night time exposure to exclude the day observations. 

We have plotted the spectra in the observed frame so that airglow features may be 
more easily identified. The intrinsic Lyman continuum region for each galaxy is also clearly 
marked. Several of the spectra are very noisy and systematic flux offsets caused by back- 
ground over-subtraction in the CALFUSE pipeline are observed. As would be expected, 
the background subtraction problems are much more prevalent in the lower S/N spectra , 
specifically SBS 0335-052, Tol 0440-381, IRAS 19245-4140, and NGC7673. A large number 
of apparent emission features can be observed, particularly in the wavelength region between 
920-930 A. These are easily observed in the day/night spectrum of VV114 but also in the 
night only spectra of Tol 0440-381. These features, while varying widely in relative strength, 
are found at the same wavelengths in most of the spectra. They are properly identified as 
geocoronal emission features, mostly H I, O I, N I, and N2. For comparison, LWRS spectra 
(SiC IB) of airglow and the earth limb are found in Figure HHl The large number of variable 
emission features highlights the difficulty in studying Lyman continuum emission in this 
wavelength regime. 

No convincing evidence of Lyman continuum emission is seen in any of the nine spectra. 
All of the galaxies have flat spectra in the Lyman continuum regions with continuum levels 
at or below zero. The spectra with continuum levels below zero are indicative of background 
over-subtraction by the CALFUSE pipeline. This is commonly observed in low S/N FUSE 
spectra on the SiC IB and SiC 2A detectors as discussed previously, and highlights the 
difficulty in studying Lyman continuum emission. The FUSE background dominates the 
signal in these wavelength regions. The emission features that are observed can all be 
identified as geo-coronal emission, with varying strengths in all of the spectra. In addition to 
these extracted spectra, we have also examined the 2-D detector images. A visual inspection 
of the 2-D spectra confirm our analysis of the one-dimensional spectra as we do not observe 
any evidence of Lyman continuum emission. 

We have chosen to derive numerical limits for the five galaxies with S/N>5 on the LiF lA 
channel. We have excluded the lower S/N galaxies such as Tol 0440-381 and IRAS 19245-4140 
due to the significant errors introduced by the background subtraction. IRAS 19245-4140 for 
example, has a S/N of only 1.3 longward of the Lyman limit on the SiC 2A channel. While 
these errors exist even in the five galaxies with higher S/N, we can still place meaningful 
limits on the Lyman escape fraction. 



We follow the work of IShapley et al.l (120061 ) in estimating both the absolute and relative 
escape fractions. The absolute escape fraction describes the fraction of ionizing photons that 
escape the host galaxy. The relative escape fraction neglects the effects of dust extinction 
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and is defined as the perc entage of escaping io nizing; pliotons relative to escaping non-ionizing 
FUV pliotons (1500 A, IShapley et al.ll2006l ). The relative escape fraction was introduced 
due to the difficulty in determining the intrinsic ionizing flux for high-z galaxies. 

For these five galaxies we have measured an upper limit to the flux below the Lyman 
limit in each galaxy. We h ave then converted these to luminosities with Lcsc ~ AFa (A=900 
A). From Starburst99 v5.1 ( Vazquez fc Leithereiil2005 ) we find an intrinsic ratio of the bolo- 
metric to ionizing radiation of Lboi/Ljoniz ~ 7. If we assume that the IR + UV luminosities 
dominate the bolometric luminosity as expected in these star-forming galaxies we derive: 



f, 



esc.abs 



Lesc / Lbol 

LuV + LiR V Lioniz 



(3) 



We slightly modify fesc.rei of IShapley et al.l (120061 ) to use 1150 A instead of 1500 A as 
our reference wavelength. These values should be roughly comparable as we would expect 
the intrinsic luminos ities of starburst galaxies to be comparable at 1150 A and 1500 A 
(ILeitherer et al.ll2002l ). We estimate fesc,rei from 
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We obtain an intrinsic ratio L ° /L » ~ 4 from StarburstQD. The derived values 

1150A' 900A 

for fesc rei and fesc abs are found in T able [231 The observed fluxes and luminosities have been 



corrected for MW extinction using ICardelli et al.l (119891 ). We find fesc.abs < 0.2 — 7% and 
fesc,rei < 9 — 26% for the galaxies. The weakest limit, fp sn.rei < 26%, is from y V114 which 
is known to contain strong intrinsic dust attenuation (iGoldader et al.ll2002l ). Our upper 
limits for fpsr,rp i are c omparable to the average value of 14 LBGs with < fesc.rei >~ 14% from 
Shapley et al.l (120061 ). This average value in the high redshift sample however is dominated 
by two galaxies with fesc.rei > 40%. While our low redshift limits are compatible with the 
sample averaged values of fesc.rei in LBGs, our hmits show that the extreme Lyman continuum 
leakage observed in some LBGs is not observed in our small, local sample. 



3. 5. 2. Limits from C II 



Given the difficulties in directly observing Lyman continuum emission on the SiC detec- 
tors it is useful t o invest igate the Lyman escape fraction using a different method. Previously, 
Heckman et al.l (l2001bl ) have argued that since the optical depth at the Lyman edge must 
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be much larger than that of the optical depth at line center for C II, we can then use the 
residual intensity in the core of the C II line to place upper limits on the escape fraction at 
the Lyman limit. 



Assuming that C II is the dominant ionic species of carbon in the H I phase, iHeckman et al 



mom derive 

TLy = 4 X 10-1^ (Nh/Nj) (b/fX) Tc = 8.6 Z^' (6/100 km s'^) tc (5) 

where Zq is the gas-phase carbon abundance in solar units and b is the Doppler parameter. 
In our earlier gaussian fits to C II, we find that b varies from ~ 65 — 470kms~^. Using the 
smallest Doppler parameter and solar abundances we expect the optical depth at the Lyman 
edge would be a factor of ~ 5 times larger than the C II optical depth. As the majority of 
the galaxies in our sample have broader lines and sub-solar abundances, we would expect 
the Lyman limit optical depth to be enhanced by an even greater value. 



An examination of the C II feature for the galaxies in our sample (e.g. Figures [20] - [351) 
shows that it is nearly black in almost all cases. Given our estimates of the optical depth 
enhancement expected at the Lyman limit, this suggests that none of the galaxies in our 
sample have significant amounts of escaping Lyman continuum emission. 

The most likely scenario for the escape of Lyman continuum emission is the picket 
fence model. Instead of a uniform slab of hydrogen surrounding a galaxy, holes in the gas 
surrounding the galaxy would allow ionizing photons to escape into the IGM. Starburst 
driven outflows could, for example, create "tunnels" in the surrounding neutral gas. The 
picket fence model would explain the variability seen in the center of some of the C II features. 
For example, in Figure [21] of VV114's C II profile a small fraction of photons are escaping 
at kms~^ while the line is completely black at both ~ ±100 km s^^. The few galaxies with 
a significant C II escape fraction stand out due to their low metallicity, suggesting that their 
Lyman continuum escape fraction would still be very low. 

We have conservatively estimated the continuum level and calculated a residual relative 
intensity for the C II line center. We calculate an upper limit to fesc,cii by finding the upper 
bound to the normalized residual flux near line center. Only our highest redshift galaxy Mrk 
54 has been excluded, as the C II feature falls in the gap between the LiF channels. We also 
estimate upper limits for Lyman continuum emission using fesc,abs= fesc.cii Luv/(Luv + Lir) 
following the methods of the previous section. For all but one of the galaxies in our sample 
we find fesc,abs< 2%. 

Using Equation [5] we have also estimated the optical depth of the Lyman limit using our 
estimate of the C II escape fraction and assuming the slab model of hydrogen surrounding 
a galaxy. The Doppler parameters are found from our Gaussian fits in Tables [3] - [TH As 
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many of these lines are non-gaussian these must be considered only rough estimates. We 
have used the gas phase oxygen metallic i ties in Table [1] for the carbon metallicity and solar 
oxygen abundances from lAsplund et al.l (120051 ) . The results can be found in Table [2^ Our 
lowest derived optical depth is ~ 20 found for NGC3310 but the majority of galaxies 
have TLy > 100. While these optical depths are only rough estimates, they suggest that the 
Lyman limit optical depth is well above the threshold needed to allow Lyman photons to 
escape the galaxy in the case of the simple slab model. 



4. Conclusions 

Starburst driven outflows are a cosmologically important example of galactic feedback. 
These outflows drive gas, metals, and energy out of galaxy centers and possibly into the IGM. 
Outflows in the early universe may also blow out "tunnels" allowing Lyman continuum pho- 
tons from stars to escape the galaxy and reionize the universe. Studies of these superwinds 
are therefore vital to expanding our understanding of galaxy formation and evolution. 

The extreme scale and energetics of starburst driven outflows complicate studies of 
their properties. Originating as a hot (T > 10*^ K) wind fluid, the outflows expand out into 
ambient gas of the galaxy and possibly beyond. Shocks, other hydrodynamical interactions, 
and entrained ambient gas and dust produce a complicated, multiphase outflow that requires 
panchromatic observations to properly characterize it. 

Observations in the FUV play an important role in understanding outflows as they 
probe the neutral, ionized, and coronal gas. This gas traces the interaction between the 
hot outflowing wind fluid and the denser ambient ISM. Unfortunately, due to observational 
constraints at both low redshifts (due to UV atmospheric absorption) and high redshifts (due 
to contamination by the Lya forest), studies at these wavelengths were severely limited until 
the launch of the FUSE satellite in 1999. As the number of star-forming galaxies observed 
by FUSE has grown, it now possible to study the global properties of these outflows in the 
FUV. 

In this work we have analyzed FUV spectra of a sample of 16 local, star-forming galax- 
ies covering almost three orders of magnitude in SFR and over two orders of magnitude 
in stellar mass. Interstellar medium absorption and stellar photospheric features are ob- 
served throughout the sample of spectra. The strongest ISM absorption lines are generally 
blueshifted by ~ 80 to 300 kms~^ relative to thge galaxy systemic velocity, implying that 
galactic outflows are the norm in starbursts. 

Coronal gas is an important coolant for gas at temperatures T~ 10^'^ K and is only 
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directly observable in the FUV through observations of the O VIAA1032, 1038 doublet. O VI 
emission is only detected in 1 of the 16 galaxies in our sample. This suggests that radiative 
cooling of the coronal gas is not dynamically significant in the wind and does not impede 
its outflow. In contrast, O VI absorption is nearly ubiquitous in our sample and is found to 
trace systematically higher outflow speeds than the photoionized and neutral gas. We find a 
strong relationship between O VI column density and the width of the line that is consistent 
with the production of O VI in a radiatively cooling gas flow passing through the coronal 
temperature regime. This is expected to occur at the interface between the outrushing hot 
wind seen in X-rays and the cool dense material in its path. 



Previous studies (I Vazquez et al.l l2004l : ISchwartz fc Martini l2004l : iGrimes et al.l |2006| ) 

have suggested that the observed outflow velocities in a galaxy increase with the ioniza- 
tion potential of the species observed. While this is clearly the case for the coronal phase 
gas, it is not the case for the neutral and photoionized material. Our new results show 
that the derived outflow velocities in thius gas in a given galaxy increase with line strength 
(equivalent width) and not ionization state. This could be explained if the observations of 
the stronger lines probe lower column density gas than the weaker absorption features. The 
lower column density gas extends to greater outflow velocities and is not readily observable 
in weaker absorption features. We conclude that the extensive studies of outflows in the 
local star forming galaxies using the Na IAA5890, 5896 doublet in the optical should be able 
to reliably measure the outflow velocities that would be seen in the strong FUV absorption- 
lines. This is because the Na IAA5890, 5896 hne should be sufficiently strong to probe the 
low column density outflowing gas. 



In agreement with previous U V (ISchwartz et al.ll2006l ) and optical (jMartinll2005l : iRupke et al 



20051 ) studies we flnd that the outflow velocities of the strongest absorption lines increase 
with increasing SFR and increasing speciflc SFR (SFR/M*). We speculate that starbursts 
with higher SFR generate more wind momentum which can accelerate the absorbing mate- 
rial (clouds) entrained in the wind to higher outflow velocities. Our sample does not probe 
the regi me of very h igh SFR (of-order 10^ M© per year) where the relationship appears to 
flatten jMartinlboOsh . 



We find that the equivalent width of the optically thin lines arising in the HI phase 
correlates most strongly with the gas-phase metallicity (as expected since the equivalent 
width should be proportional to the column density of the metal species). In contrast, the 
equivalent widths of the optically-thick lines in the HI and HII phase correlate better with 
the galaxy mass and the SFR. This is consistent with the fact that the equivalent width for 
these optically-thick lines is determined by the velocity dispersion in the ISM. More massive 
galaxies and galaxies undergoing elevated SFR will have higher velocity dispersions due to 
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both gravity and the kinetic energy input from massive stars. 

Lastly, we have looked for escaping Lyman continuum emission in our sample. We find 
no convincing evidence of escaping ionizing radiation in any of the observed galaxies. In 
a subsample of five of our galaxies with suitable redshifts and data quality we find that 
the upper limit to the absolu te fraction of e s capin g Lyman continuum photons ranges from 



fesc,abs< 0.2 — 6.3%. Following lShaplev et al.l (120061 ) we also derive upper limits to the realtive 



escape fractioni (the ratio of the fraction of escaping Lyman continuum and non-ionizing FUV 
continuum). These range from fesc,rei< 9 — 26% for this subsample. While our upper limits 
are compatible with observations of high redshift Lyman Break Galaxies, we do not observe 
any evidence of Lyman continuum emission as observed in some Lyman Break Galaxies. 

Using the escape fraction of light at the center of the C II feature we also place limits 
on the escape fraction of Lyman continuum emission in 15 of the galaxies included in our 
sample. In a simple model where a galaxy is surrounded by a slab of neutral hydrogen we 
find significant optical depths at the Lyman limit of > 20. In a more realistic model 
with holes in the neutral gas surrounding a galaxy, we find fesc,abs< 2% for almost all of the 
galaxies. This suggests that local starburst galaxies are not allowing the escape of significant 
amounts of hydrogen ionizing photons, even in galaxies with strong winds. 

This work provides a unique view of the properties of outflows in star-forming galaxies. 
The wealth of absorption features over a variety of gas phases enables direct comparisons dif- 
ficult in other wavelength regimes. The projected physical size of FUSE^s large spectroscopic 
aperture for these local starbursts is particularly well matched to high redshift spectrographic 
observations of Lyman Break Galaxies and other star forming galaxies - both probe the out- 
flows on the galactic scale. Thus, this study should provide an excellent complement to the 
high redshift observations. 
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IRAS19245-4140 


2832 


40 


0.087 


8.8 


X 


109 


3.9 


X 


10« 


8.2 


X 


109 


1.7 


X IQio 


1.67 


8.0 


NGC7673 


3392 


48 


0.043 


6.7 


X 


109 


3.2 


X 


109 


3.2 


X 


IQio 


3.9 


X IQio 


0.74 


8.5 


NGC7714 


2803 


39 


0.052 


9.4 


X 


108 


5.3 


X 


109 


5.6 


X 


IQio 


5.7 


X IQio 


0.91 


8.7 



00 
00 



No t e. — Vsva, MW E(B -V), IRAS, and K-Band flux e s wer e obtained from NED where available. IRAS jSoifer. ' joehmgrj^jeugeba^ Sanders! 
19891: iMoshir et al.l Il993h and 2Mass lljarrett et al.l l2003l) luminosities are transformed using methods of ISande^^^Mirab cl ' /l99 ^ and 
Carpente A l|200l|}_ respectively. The references for the gas -phase nebular metallicity measurements are as follows, Haro 11 ijBcrgva ll fc 6stlinl 



200 j). VV 114 jKim et "Zl Il995l : Icharlot fc Longhettil [2OO1I) . NGC 1 140, NGC1705, NGC 3310, NGC 3690, NG C4214, NG C 7673, NGC 7714 



llHeckmanetal, 
Haro 3 lllzotov ~ 



- - - gti , - , , 

1998h. SBS 335-052 lllzotov et al.lll997[). Tol 0440-381 iFerlevich et al.lll99ll: jPeilerin fc Rober3l2007h . NGC 1741 dConti et al.lll996l) . 
Thuanlliooi). IZwl8|Izotov fc Thuanlll999l) . IRAS 19245-4140 llOstlin et allbOOTi) . and Mrk54 (spectrum from A. Boselh with 



versions from lcharioT fc LonghettilboOlir 

'^FU V luminosity derived from FUSE LiF 2A continuum flux at 1150 A. Values are corrected for MW extinction using the method of lCardelli et al.l 
lll989h and the E(B-V) values given in this table. 



Table 2. FUSE Observation Data 



Galaxy 


FUSE Pointing 
(J2000) 


Aperture 


• 

Region 




L/Ds. JJate 




\Jos. lime 
(ksj 


Exposures 


Q /MC 
o/lN 


HARO 11 


00''36™52?5 


-33°33'19" 


LWRS 


PC 


NIGHT 


2001-10-12 


Bergvall 


11.9 


6 


7.2 


VV 114 


0lh07m4gsg 


-17°30'24" 


LWRS 


PC 


BOTH 


2003-07-26 


Hcckman 


11.2 


4 


6.6 


NGC 1140 


02''54™33?6 


-10°01'39" 


LWRS 


PC 


BOTH 


2001-11-28 


Buat 


3.9 


2 


4.8 


SBS 0335-052 


03'^37™43!1 


-05°02'39" 


LWRS 


PC 


BOTH 


2001-09-26 


Thuari 


24.4 


8 


4.1 


TOL 0440-381 


04'i42™08?0 


-38°01'10" 


LWRS 


PC 


BOTH 


2000-12-13 


Deharveng 


31.2 


11 


4.0 


NGC 1705 


04''54™13?5 


-53°21'39" 


LWRS 


EC 


BOTH 


2000-02-04 


Heckman 


24.0 


13 


27.2 


NGC 1741 


05''01™37!7 


-04° 15' 30" 


LWRS 


PC 


NIGHT 


2003-12-23 


Heckman 


18.5 


10 


8.4 


IZW 18 


09''34™02?3 


-|-55°14'25" 


LWRS 


PC 


BOTH 


2001-02-11 


Moos 


58.2 


37 


9.7 


NGC 3310 


10'^38™45?7 


-l-53°30'06" 


LWRS 


EC 


BOTH 


2000-05-05 


Heckman 


27.1 


11 


18.0 


HARO 3 


10'^45™22?4 


-l-55°57'37" 


LWRS 


PC 


NIGHT 


2002-04-07 


Thuan 


4.9 


4 


6.2 


NGC 3690 


ll'^28™31?0 


-l-58°33'41" 


LWRS 


PC 


NIGHT 


2001-03-25 


Heckman 


34.9 


20 


10.8 


NGC 4214 


12''15™39?4 


-1-36° 19' 35" 


LWRS 


EC 


BOTH 


2000-05-12 


Heckman 


11.1 


6 


14.6 


MRK 54 


12''56™55!9 


-|-32°26'52" 


LWRS 


PC 


BOTH 


2000-02-19 


Deharveng 


25.8 


13 


5.8 


IRAS 19245-4140 


19'^27™58!0 


-41°34'27" 


LWRS 


PC 


NIGHT 


2000-05-21 


Buat 


2.9 


6 


4.7 


NGC 7673 


23'^27™41!6 


-|-23°35'30" 


LWRS 


PC 


NIGHT 


2000-07-20 


Buat 


6.3 


3 


4.7 


NGC 7714 


23'' 36™ 14! 1 


-|-02°09'19" 


MORS 


PC 


BOTH 


2001-12-04 


Benjamin 


78.4 


46 


4.3 



^Detector region used for spectra extraction. Although all targets are extended (EC) the majority have flux falling entirely within the smaller point 
source (PC) extraction region. Using the smaller point source region minimizes background subtraction errors which are very important in the lower S /N 
data. 

''Exposure time in LiF lA channel. 
<=S/N in 0.078 A bins in LiF lA channel. 
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Table 3. HAROll ISM Absorption Line Fit Data 



Ion 




loglXfN/NH) 


Instrument 




Vc 


FWHM 




A 






A 


kms~' 


kms""' 


cm 


977.02 


-0.61 


2asic 


1.83±0.57 


5998±6 


269±10 


oVI 


1031.93 


-1.13 


lalif 


0.38±0.06 


6039±16 


193±34 








2blif 


0.39±0.07 


6050±14 


159±34 


TT 

C 11 


1036.34 


-1.37 


lalif 


1.29±0.05 


6063±6 


381±18 








2blif 


1.24±0.13 


6074±8 


344±36 


ol 


988.73 


-1.56 


lalif 


0.72±0.05 


6092±14 


276±14 








2blif 


0.66±0.09 


6079±23 


304±26 


Nil 


1083.99 


-2.00 


Iblif 


0.92±0.07 


6081±10 


353±28 








2alif 


1.00±0.08 


6085±11 


381±33 


NlII 


989.80 


-2.03 


lalif 


1.04±0.05 


6076±6 


276'' 








2blif 


1.05±0.08 


6073±10 


304'" 


Ol 


1039.23 


-2.29 


lalif 


0.35±0.07 


6112±14 


244±51 








2blif 


0.37±0.12 


6130±16 


200±77 


Nl 


1134.41 


-2.57 


2alif 


Cl-Si" 


6048±79 


300'' 


Fell 


1096.88 


-2.96 


2alif 


0.20±0.04 


6119±9 


131±24 


Si II 


1020.70 


-2.99 


lalif 


0.53±0.18 


6023±17 


339±88 








2blif 


0.32±0.10 


6015±17 


232±67 


Fell 


1121.97 


-3.14 


2alif 


0.25'' 


6088±68 


203'' 


SlII 


1012.50 


-3.24 


lalif 


0.52±0.18 


6083±6 


225'' 








2blif 


0.68±0.22 


6078±13 


290±75 


Fell 


1125.45 


-3.41 


2alif 


0.23±0.05 


6062±17 


195±37 


Arl 


1066.66 


-3.58 


2alif 


0.54±0.08 


6136±20 


332±55 




1128.00 


-4.01 


2alif 


0.14*= 


6202±74 


327'' 


pVs 


1117.98 




2alif 


0.30±0.06 


6165±24 


327±63 


SilVs 


1122.48 




2alif 


0.63±0.09 


6163±17 


222±29 


SilVs 


1128.33 




2alif 


0.30'' 


6202'' 


327'' 



''Using Vsys = 6175 kms ' 

''Errors are undetermined for this value. 
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Table 4. HAR03 ISM Absorption Line Fit Data 



Ion 


Ao 




Instrument 






FWHM 




A 






A 


km s ~ 


km s^"*^ 


ClII 


977.02 


-0.61 


2asic 


1.16*^ 


913±26 


202±16 


oVI 


1031.93 


-1.13 


lalif 


0.71±0.08 


787±26 


324'' 








2blif 


0.85±0.12 


796±21 


324±53 


cll 


1036.34 


-1.37 


lalif 


0.97±0.07 


868±3 


229±26 








2blif 


1.11±0.13 


876±10 


259±54 


ol 


988.73 


-1.56 


lalif 




912±11 


150'' 








2blif 


0.46±0.06 


906±4 


88±14 


NlII 


989.80 


-2.03 


lalif 


0.87'' 


914±6 


127it2 








2blif 


0.59'' 


920±22 


no'' 


ol 


1039.23 


-2.29 


lalif 


0.51±0.04 


895±4 


120±16 








2blif 


0.47±0.10 


905±2 


86±5 


Nl 


1134.98 


-2.39 


Iblif 


0.33±0.04 


900±8 


118'' 








2alif 


0.30±0.04 


901±12 


144b 


Fell 


1144.94 


-2.42 


Iblif 


0.42±0.07 


891±5 


85±22 








2alif 


0.45±0.10 


871±12 


172it44 


Nl 


1134.41 


-2.57 


Iblif 


0.41±0.04 


884±6 


118±11 








2alif 


0.38±0.04 


878±8 


144±14 


Fell 


1096.88 


-2.96 


Iblif 


0.34±0.06 


904±8 


117±20 








2alif 


0.31±0.06 


886dz9 


139±24 


Sill 


1020.70 


-2.99 


lalif 


0.35±0.04 


912±4 


107±15 








2blif 


0.41±0.13 


941±15 


118±37 


Fell 


1121.97 


-3.14 


Iblif 


0.43±0.08 


902±9 


132±29 








2alif 


0.46±0.12 


924±16 


156±42 


Fell 


1143.23 


-3.19 


Iblif 


0.21±0.08 


894±12 


86'' 








2alif 


0.25±0.07 


907±13 


134±36 


SlII 


1012.50 


-3.24 


lalif 


0.40±0.07 


921±8 


130±22 








2blif 


0.45it0.06 


920±6 


109±18 


Arl 


1066.66 


-3.58 


lalif 


0.14'' 


888±31 


80'' 








2blif 


0.19±0.07 


897±10 


67'' 


Si IV s 


1122.48 




Iblif 


0.47±0.08 


922±7 


129±26 








2alif 


0.44±0.06 


902±4 


105±15 


ClIIs 


1175.59 




Iblif 


1.50it0.15 


938±31 


500'' 



'■Using Vsys = 963 kms ^ 

''Errors are undetermined for this value. 



-37- 



Table 5. IRAS19245-4140 ISM Absorption Line Fit Data 



Ion 


Ao 

A 


log{\fN/NH) 


Instrument 


A 


Vc 


FWHM 

kms^-' 


cm 


977.02 


-0.61 


2asic 


2.50±0.74 


2748±29 


663it235 


ClI 


1036.34 


-1.37 


lalif 


0.86±0.09 


2746±11 


272±30 








2blif 


0.75±0.25 


2761±31 


256±96 


ol 


988.73 


-1.56 


lalif 


0.31±0.03 


2774±8 


81±13 


Nil 


1083.99 


-2.00 


2alif 


0.17^ 


2792±14 


72'' 


NlII 


989.80 


-2.03 


lalif 


0.32^ 


2708±7 


89'' 


Nl 


1134.41 


-2.57 


Iblif 


0.38*^ 


2782±1 


50±11 








2alif 


0.30*= 


2794±0 


37'' 


Si II 


1020.70 


-2.99 


lalif 


0.14±0.05 


2821±9 


59±21 








2blif 


0.07'' 


2820±27 


59'' 


Arl 


1048.22 


-3.00 


lalif 


0.08'' 


2787±14 


48'' 


Fell 


1121.97 


-3.14 


Iblif 


O.ll'' 


2862±15 


67'' 








2alif 


0.14'' 


2828±11 


45'' 


SlII 


1012.50 


-3.24 


lalif 


0.32^ 


2785±5 


50±7 








2blif 


0.30*' 


2793±8 


50'' 


SilVs 


1122.48 




Iblif 


0.18±0.02 


2830±6 


60±12 








2alif 


0.21'' 


2829±11 


60'' 


ClIIs 


1175.59 




Iblif 


2.55±0.50 


2795±25 


641±130 



aUsing vsys = 2832 kms"! 

''Errors are undetermined for this value. 
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Table 6. IZW18 ISM Absorption Line Fit Data 



Ion 


Ao 


log{\fN/NH) 


Instrument 




Vc 


FWHM 




A 






A 






ClI 


1036.34 


-1.37 


lalif 


0.63±0.05 


748±7 


180±16 








2blif 


0.65±0.06 


757±8 


190±18 


Ol 


988.73 


-1.56 


lalif 


0.64±0.23 


784±4 


106±12 








2asic 


0.54'> 


816±5 


84±15 


NlII 


989.80 


-2.03 


lalif 


0.51±0.10 


711±9 


122±28 








2asic 


0.69±0.15 


754±11 


180±68 


Ol 


1039.23 


-2.29 


lalif 


0.18±0.02 


759±5 


76±9 








2blif 


0.23±0.03 


770±5 


79±13 


Si II 


1020.70 


-2.99 


lalif 


0.15±0.02 


774±8 


lOO'' 








2blif 


0.17±0.03 


788±11 


lOO'' 


Fell 


1121.97 


-3.14 


Iblif 


0.10±0.03 


751±8 


65±22 








2alif 


0.13±0.01 


769±2 


57±6 


Fell 


1125.45 


-3.41 


Iblif 


0.09±0.02 


745±6 


56±16 








2alif 


0.11±0.02 


772±6 


79±15 


Arl 


1066.66 


-3.58 


lalif 


0.06'' 


758±20 


70'' 


Si IV s 


1122.48 




Iblif 


0.18±0.05 


753±17 


149±54 








2alif 


0.21±0.02 


768±5 


126±15 


SilVs 


1128.33 




Iblif 


O.IS'' 


775±34 


155'' 








2alif 


O.IS*" 


779±68 


228*' 


ClIIs 


1175.59 




Iblif 


0.86±0.07 


747±15 


423±38 








2alif 


0.81±0.06 


757±12 


388±30 



■^Using Vays = 751 km s ^ 

''Errors are undetermined for this value. 
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Table 7. MRK54 ISM Absorption Line Fit Data 



Ion 


Ao 


log{\fN/NH) 


Instrument 




Vc 


FWHM 




A 






A 






cm 


977.02 


-0.61 


lalif 


1.76±0.12 


13382±12 


434it46 








2blif 


1.48±0.09 


13402±6 


368±28 


oVI 


1031.93 


-1.13 


lalif 


0.77±0.09 


13321±14 


287±38 


ol 


988.73 


-1.56 


lalif 


1.18±0.34 


13541±46 


386±112 








2blif 


0.82±0.19 


13526±30 


337±86 


Nil 


1083.99 


-2.00 


Iblif 


1.41±0.12 


13466±9 


409±39 








2alif 


1.60±0.13 


13469±10 


497±44 


NlII 


989.80 


-2.03 


lalif 


0.67±0.11 


13481±13 


177±45 








2blif 


0.82'' 


13463±5 


138±20 


Si II 


1020.70 


-2.99 


lalif 


0.53±0.15 


13544±58 


300'' 








2blif 


0.81±0.07 


13528±21 


400*' 


SlII 


1012.50 


-3.24 


lalif 


0.62±0.11 


13432±14 


310±46 








2blif 


0.54±0.12 


13429±15 


277±56 


SiVs 


1062.66 




Iblif 


0.64±0.07 


13529±27 


400'' 








2alif 


0.74±0.06 


13525±20 


400'' 


SilVs 


1122.48 




Iblif 


0.96±0.12 


13434±14 


312±44 








2alif 


0.97±0.16 


13423±24 


334±62 


SilVs 


1128.33 




2alif 


0.71±0.07 


13453±14 


299±33 



''Using Vsys = 13450 kms~' 

''Errors are undetermined for this value. 
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Table 8. NGC1140 ISM Absorption Line Fit Data 



Ion 




/ogfAf .V/.Vrr) 


Iiistruiiiciit 


Wa 


Vc 


FWHM 




A 






A 


kms 1 


kms-l 


ClII 


977.02 


-0.61 


2asic 


1.32±0.18 


1438±11 


268±46 


ClI 


1036.34 


-1.37 


lalif 


0.93±0.08 


1481±6 


250±21 








2blif 


1.33±0.27 


1487±27 


358±76 


ol 


988.73 


-1.56 


2asic 


0.51±0.12 


1542±5 


100'' 








2blif 


O.TS'^ 


1548±6 


lOO'' 


Nil 


1083.99 


-2.00 


2alif 


0.63±0.10 


1517±7 


177±28 


NlII 


989.80 


-2.03 


2asic 


0.99'' 


1502±5 


138±39 








2blif 


0.69±0.10 


1517±7 


130'' 


Ol 


1039.23 


-2.29 


lalif 


0.45±0.03 


1510±2 


llOill 








2blif 


0.42±0.05 


1526±4 


112±13 


Nl 


1134.98 


-2.39 


Iblif 


0.24±0.05 


1505±12 


109±21 








2alif 


0.22±0.03 


1515±8 


103'' 


Fell 


1144.94 


-2.42 


Iblif 


0.39±0.04 


1511±6 


105'' 








2alif 


0.39±0.06 


1534±6 


121±21 


Nl 


1134.41 


-2.57 


Iblif 


0.30±0.04 


1490±8 


109±14 








2alif 


0.26±0.03 


1491±7 


103±11 


Fell 


1063.18 


-2.70 


lalif 


0.37±0.06 


1504±9 


145±30 








2blif 


0.44±0.09 


1514±13 


147±37 


Fell 


1096.88 


-2.96 


Iblif 


0.26±0.05 


1503±9 


122±28 








2alif 


0.27±0.04 


1515±7 


112±21 


Arl 


1048.22 


-3.00 


lalif 


0.16±0.04 


1512±16 


130'' 








2blif 


0.26±0.05 


1529±13 


130'' 


Fell 


1121.97 


-3.14 


Iblif 


0.31±0.10 


1512±10 


115±42 








2alif 


0.32'' 


1494±25 


148'' 


SlII 


1012.50 


-3.24 


lalif 


0.47±0.09 


1509±9 


175±31 








2blif 


0.50'' 


1501±27 


168'' 


Arl 


1066.66 


-3.58 


lalif 


0.16±0.04 


1493±24 


130'' 


SilVs 


1122.48 




Iblif 


0.41±0.08 


1512±10 


132±31 








2alif 


0.67±0.18 


1494±25 


208±47 


ClIIs 


1175.59 




Iblif 


1.16±0.07 


1498±14 


300'' 



''Using Vsys = 1501 km s ' 

''Errors are undetermined for this value. 
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Table 9. NGC1705 ISM Absorption Line Fit Data 



Ion 


Ao 


log{X}N/NH) 


Instrument 




Vc 


FWHM 




A 






A 


kms~' 


kms~' 


/-I TTT 
C ill 


977.02 


-0.61 


Ibsic 


1.21it0.29 


613±20 


310 








2asic 


1.23'= 


651it28 


325'= 


O V i 


1031.93 


-1.13 


lalif 


0.18±U.01 


554i:3 


lllit9 








2blif 


O.iOiO.Ol 


559±3 


125±9 


<-i TT 


1036.34 


-1.37 


lalif 


1.01*= 


600±37 


262'' 








2blif 


1.02±0.2d 


597±24 


270±84 


r\ T 
O 1 


988.73 


-1.56 


lalif 


0.58±0.04 


600±3 


162±21 








2blif 


0.53±0.04 


596±4 


153±14 


M TT 
N 11 


1083.99 


-2.00 


lasic 


0.49±0.0d 


626±4 


160±18 


AT TTT 

N ill 


989.80 


-2.03 


lalif 


0.60±0.05 


597±2 


130±13 








2blif 


0.63±0.04 


594±3 


145±14 


r\ T 
O 1 


1039.23 


-2.29 


lalif 


0.22*= 


606±31 


72'= 








2blif 


0.23'' 


600±32 


77'' 


Nl 


1134.98 


-2.39 


Iblif 


0.07±0.01 


592±6 


73'= 








2alif 


0.05±0.01 


593±3 


50'= 


Fell 


1144.94 


-2.42 


Iblif 


0.24±0.02 


610±2 


87±7 








2alif 


0.25±0.02 


609±2 


89±7 


Nl 


1134.41 


-2.57 


Iblif 


0.07±0.01 


614±9 


73±9 








2alif 


0.05±0.01 


597±8 


50±15 


Fell 


1096.88 


-2.96 


Iblif 


0.14±0.01 


607±3 


88±7 








2alif 


0.11±0.02 


603±5 


67±14 


Si II 


1020.70 


-2.99 


lalif 


0.09±0.01 


613±2 


77±7 








2blif 


0.09±0.01 


599±3 


80'= 


Fell 


1121.97 


-3.14 


Iblif 


0.18'= 


620±50 


120'= 








2alif 


0.12'' 


603±7 


108±43 


SlII 


1012.50 


-3.24 


lalif 


0.18±0.02 


608±11 


130'= 








2blif 


0.17±0.03 


603±14 


130'= 


Arl 


1066.66 


-3.58 


lalif 


0.11±0.02 


622±9 


110'= 








2blif 


0.12±0.01 


592±6 


110'= 


SilVs 


1122.48 




Iblif 


0.40'' 


627±47 


172'= 








2alif 


0.48±0.00 


618±4 


200±4 


ClIIs 


1175.59 




Iblif 


1.67±0.04 


639±4 


515±13 








2alif 


1.65±0.04 


634±4 


507±12 



'■Using Vgys = 632 kms ^ 

''Errors are undetermined for this value. 



-42 - 



Table 10. NGC1741 ISM Absorption Line Fit Data 



Ion 


Ao 


log{\fN/NH) 


Instrument 




Vc 


FWHM 




A 






A 






ri TTT 
C ill 


977.02 


-0.61 


2asic 


1.08 


4007±20 


164±15 


O Vi 


1031.93 


-1.13 


lalif 


0.58±0.08 


3833it32 


400'' 








2blif 


0.58±0.20 


3921±92 


400*' 




1036.34 


-1.37 


lalif 


1.38±0.07 


3985±4 


300±28 








2blif 


1.39±0.19 


3986±4 


270±29 


Ol 


988.73 


-1.56 


lalif 


1.16±0.25 


4024±10 


234±56 








2blif 


0.98*^ 


4046±11 


182±50 


Nil 


1083.99 


-2.00 


Iblif 


0.93±0.19 


3957±9 


277±57 








2alif 


0.91±0.05 


3982±5 


263±16 


NlII 


989.80 


-2.03 


lalif 


0.84±0.10 


4004±12 


223±31 








2blif 


0.79±0.15 


3981±18 


206±37 


Ol 


1039.23 


-2.29 


lalif 


0.58±0.07 


3981±7 


222±25 








2blif 


0.49±0.07 


3977±11 


190±28 


Nl 


1134.98 


-2.39 


Iblif 


0.16±0.06 


3962±11 


82±24 








2alif 


0.16±0.04 


3944±16 


141±40 


Fell 


1144.94 


-2.42 


2alif 


0.40±0.04 


3992±13 


200'' 


Nl 


1134.41 


-2.57 


Iblif 


0.20±0.07 


3937±9 


82±31 








2alif 


0.31±0.05 


3928±15 


141±23 


Fell 


1096.88 


-2.96 


Iblif 


0.30±0.05 


3953±13 


181±37 








2alif 


0.30±0.06 


3964±13 


200±46 


Si II 


1020.70 


-2.99 


lalif 


0.45±0.08 


3954±8 


211±34 








2blif 


0.42±0.08 


3963±8 


180±35 


Fell 


1121.97 


-3.14 


Iblif 


0.45±0.14 


3979±22 


204±54 








2alif 


0.34^ 


3978±42 


199±70 


SiVs 


1062.66 




lalif 


0.38±0.05 


4017±23 


200'' 


Si IV s 


1122.48 




Iblif 


0.58±0.10 


4001±13 


194±29 








2alif 


0.59±0.15 


4003±23 


218±39 


''Using 


vsys = 4037 kms-l 











'Errors are undetermined for this value. 
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Table 11. NGC3310 ISM Absorption Line Fit Data 



Ion 




/ogfAf .V/.Vrr) 


Iristriirricnt 




Vc 


FWHM 




A 






A 


kms ^ 


km s - 1 


ClII 


977.02 


-0.61 


Ibsic 


2.44±0.05 


701±5 


583±19 








2asic 


2.10±0.61 


705±19 


450±169 


oVI 


1031.93 


-1.13 


lalif 


1.67±0.04 


616±5 


556±15 








2blif 


1.71±0.04 


620±5 


554±15 


ClI 


1036.34 


-1.37 


lalif 


3.31±0.07 


749±6 


786±31 








2blif 


3.41±0.07 


752±6 


803±35 


ol 


988.73 


-1.56 


lalif 


1.57'= 


877±10 


270±50 








2blif 


1.52'' 


859±6 


259±42 


Nil 


1083.99 


-2.00 


lasic 


2.23±0.17 


889±9 


627±48 








2bsic 


2.06±0.13 


874±9 


581±40 


NlII 


989.80 


-2.03 


lalif 


1.21±0.13 


891±17 


349±38 








2blif 


1.11±0.15 


884±21 


352±47 


Ol 


1039.23 


-2.29 


lalif 


0.63±0.05 


945±16 


450'' 








2blif 


0.72±0.01 


966±15 


450'' 


Nl 


1134.98 


-2.39 


Iblif 


0.16±0.04 


837±6 


397±18 








2alif 


0.15±0.04 


846±7 


360±19 


Nl 


1134.41 


-2.57 


Iblif 


0.55±0.03 


837±6 


397±19 








2alif 


0.50±0.03 


846±7 


360±19 


Fell 


1096.88 


-2.96 


2alif 


0.40±0.11 


892±38 


351±101 


Fell 


1121.97 


-3.14 


Iblif 


0.45±0.10 


885±20 


305±48 








2alif 


0.53'' 


909±32 


335±47 


SlII 


1012.50 


-3.24 


lalif 


0.56±0.02 


950±8 


350'' 








2blif 


0.78±0.22 


923±10 


439±95 


Arl 


1066.66 


-3.58 


lalif 


0.39±0.02 


930±4 


231±15 








2blif 


0.34±0.06 


922±13 


229±39 


SiVs 


1062.66 




lalif 


0.82±0.19 


946±35 


520±106 








2blif 


0.64'' 


949±105 


432'' 


SilVs 


1122.48 




Iblif 


0.70±0.12 


978±20 


346±50 








2alif 


0.71±0.25 


973±48 


383±73 


ClIIs 


1175.59 




2alif 


1.54±0.05 


985±4 


530±15 



^Using Vsys = 983 kms~l 

''Errors are undetermined for this value. 
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Table 12. NGC3690 ISM Absorption Line Fit Data 



Ion 


Ao 


log{\fN/NH) 


Instrument 




Vc 


FWHM 




A 






A 






cm 


977.02 


-0.61 


Ibsic 


1.69±0.10 


3033±84 


385*' 








2asic 


1.53±0.20 


3111±9 


364±104 


oVI 


1031.93 


-1.13 


lalif 


0.78'' 


2866±37 


429*' 








2blif 


0.98±0.16 


2842±24 


564±89 


ClI 


1036.34 


-1.37 


lalif 


2.40±0.10 


3051±6 


580it39 








2blif 


2.20±0.08 


3052±6 


501±36 


Nil 


1083.99 


-2.00 


2alif 


2.27±0.08 


3075±5 


589±22 


Ol 


1039.23 


-2.29 


lalif 


0.63±0.07 


3169±4 


172±20 








2blif 


0.82±0.06 


3142±4 


224±17 


Nl 


1134.98 


-2.39 


Iblif 


0.42±0.04 


3153±8 


173'= 








2alif 


0.53±0.03 


3151±6 


192±10 


Nl 


1134.41 


-2.57 


Iblif 


0.55±0.04 


3132±7 


173±12 








2alif 


0.68±0.03 


3138±5 


192±10 


Fell 


1096.88 


-2.96 


Iblif 


0.49±0.08 


3131±8 


213±35 








2alif 


0.61±0.0r5 


3098±6 


233±20 


SlII 


1012.50 


-3.24 


lalif 


0.55^ 


3099±146 


319'' 








2blif 


0.57±0.10 


3114±19 


262it41 


SilVs 


1128.33 




Iblif 


0.67±0.15 


3079±27 


365±69 








2alif 


0.48^' 


3128±36 


299±88 



^Using Vsys = 3121 kms ^ 

''Errors are undetermined for this value. 
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Table 13. NGC4214 ISM Absorption Line Fit Data 



Ion 


Ao 


loo(\fN/NH) 


Instrument 


Wx 


Vc 


FWHM 




A 






A 


kms"^ 


knis~^ 


ClII 


977.02 


-0.61 


2asic 


1.04±0.33 


233±4 


194±28 


oVI 


1031.93 


-1.13 


lalif 


0.29±0.04 


212±5 


124±16 








2blif 


0.35±0.03 


225±5 


152±15 


c II 


1036.34 


-1.37 


lalif 


1.04±0.04 


273±2 


263±13 








2blif 


1.06±0.03 


279±2 


265±12 


ol 


988.73 


-1.56 


lalif 


0.47±0.14 


283±12 


110*' 








2blif 


0.50±0.15 


285±6 


110*" 


Nil 


1083.99 


-2.00 


lasic 


0.63±0.03 


289±4 


162±8 








2bsic 


0.74±0.09 


299±7 


197±28 


NlII 


989.80 


-2.03 


lalif 


0.55±0.06 


294±4 


135±24 








2blif 


0.50±0.06 


295±2 


98±12 


Ol 


1039.23 


-2.29 


lalif 


0.41±0.02 


294±2 


112±7 








2blif 


0.47±0.03 


298±2 


128±8 


Nl 


1134.98 


-2.39 


Iblif 


0.29±0.02 


305±4 


110±6 








2alif 


0.29±0.01 


301±3 


112±11 


Nl 


1134.41 


-2.57 


Iblif 


0.35±0.02 


290±3 


110±6 








2alif 


0.34±0.02 


283±3 


112±5 


Sill 


1020.70 


-2.99 


lalif 


0.28±0.02 


312±3 


104±8 








2blif 


0.25±0.02 


308±2 


93±6 


SlII 


1012.50 


-3.24 


lalif 


0.32±0.02 


293±3 


123±9 








2blif 


0.31±0.02 


293±3 


120±9 


pVs 


1117.98 




2alif 


0.27±0.05 


294±10 


290±44 


SilVs 


1128.33 




Iblif 


0.63±0.05 


298±7 


318±25 








2alif 


0.62±0.04 


304±6 


334±21 


ClIIs 


1175.59 




Iblif 


1.53±0.07 


285±9 


545±26 








2alif 


1.40±0.07 


288±6 


510±25 



''Using Vgys = 291 km s ^ 

''Errors are undetermined for this value. 
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Table 14. NGC7673 ISM Absorption Line Fit Data 



Ion 


A 


log{\fN/NH) 


Instrument 


A 


Vc 


FWHM 
kms~^ 


oVI 


1031.93 


-1.13 


lalif 


0.95±0.17 


3323±28 


386±66 








2blif 


1.07±0.15 


3323*= 


386'' 


ClI 


1036.34 


-1.37 


lalif 


2.06±0.34 


3355±30 


469±166 








2blif 


2.07±0.58 


3382±14 


350±11 


Nil 


1083.99 


-2.00 


Iblif 


1.52±0.12 


3401±15 


440±40 








2alif 


1.45±0.12 


3405±14 


411±32 


Ol 


1039.23 


-2.29 


lalif 


1.02±0.12 


3435±8 


283±29 








2blif 


0.92±0.17 


3419±10 


254±53 


Nl 


1134.41 


-2.57 


2alif 


0.47±0.06 


3404±26 


381±53 


pVs 


1117.98 




2alif 


0.45±0.05 


3417±19 


250'' 


SilVs 


1128.33 




Iblif 


0.32±0.09 


3402±33 


250'' 








2alif 


0.57±0.10 


3455±19 


259±50 



^Using vsys = 3392 kms-^ 

''Errors are undetermined for this value. 
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Table 15. NGC7714 ISM Absorption Line Fit Data 



Ion 


A 


log{\fN/NH) 


Instrument 


A 


Vc 


FWHM 

km s^-'- 


oVI 


1031.93 


-1.13 


lalif 


1.54±0.01 


2589±3 


363±6 








2blif 


1.56±0.14 


2594±16 


522±51 


ClI 


1036.34 


-1.37 


lalif 


1.94±0.35 


2739±42 


532±110 








2blif 


2.76±0.10 


2726±9 


709±30 


Nil 


1083.99 


-2.00 


2alif 


2.36±0.11 


2753±8 


628±33 


Fell 


1096.88 


-2.96 


Iblif 


0.49±0.09 


2725±11 


220±48 








2alif 


0.46±0.08 


2757±14 


349±55 


Fell 


1121.97 


-3.14 


2alif 


0.50±0.10 


2727±19 


258±36 


SlII 


1012.50 


-3.24 


lalif 


0.58'' 


2750±61 


188'' 








2blif 


0.62±0.04 


2751±9 


200'' 


Si IV s 


1122.48 




2alif 


0.80±0.08 


2766±9 


255±22 



^Using Vsys = 2803 kms"! 

''Errors are undetermined for this value. 
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Table 16. SBS0335-052 ISM Absorption Line Fit Data 



Ion 


Ao 


log{XfN/NH) 


Instrument 


Wa 


Vc 


FWHM 




A 






A 


kms"'^ 


kms~' 


oVI 


1031.93 


-1.13 


lalif 


o.ie'' 


3863±105 


155'' 








2blif 


0.13'' 


3919±82 


132'' 


ClI 


1036.34 


-1.37 


lalif 


0.47±0.08 


4040±10 


138±26 








2blif 


0.38±0.07 


4051±7 


94±19 


ol 


988.73 


-1.56 


lalif 


0.40*' 


4043±12 


76'' 








2blif 


0.39'' 


4075±10 


56'' 


Nil 


1083.99 


-2.00 


2alif 


0.23'' 


4095±45 


176'' 


NlII 


989.80 


-2.03 


lalif 


0.46±0.10 


4058±11 


115±24 








2blif 


0.51±0.18 


4093±7 


86±12 


Ol 


1039.23 


-2.29 


lalif 


0.23'' 


4066±3 


35'' 








2blif 


0.21±0.02 


4081±3 


54±3 


Nl 


1134.98 


-2. ,39 


2alif 


0.10±0.03 


'1069±4 


31'' 


Nl 


11 


-2.57 


2alil 


().il±o.o:-! 


ll)7:-!±r, 


:ii±9 



''Using vsys = 4043 kms~' 

''Errors are undetermined for this value. 
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Table 17. TOL0440-381 ISM Absorption Line Fit Data 



Ion 


Ao 


log{\fN/NH) 


Instrument 




Vc 


FWHM 




A 






A 


kms~^ 


kms~^ 


cm 


977.02 


-0.61 


lalif 


1.23±0.22 


12063±16 


369±76 








2blif 


1.23*" 


12063±7 


369*' 


ClI 


1036.34 


-1.37 


lalif 


1.25it0.42 


12171±9 


222it31 


Nil 


1083.99 


-2.00 


2alif 


0.60±0.07 


12194±11 


206±25 


SlII 


1012.50 


-3.24 


lalif 


0.52±0.05 


12197ztl2 


252±26 








2blif 


0.58±0.06 


12205±9 


192±22 


SiVs 


1062.66 




2alif 


0.55±0.11 


12254±32 


376±73 


Si IV s 


1128.33 




2alif 


0.82±0.13 


12298±28 


351±54 



*Using Vsys = 12291 kms"! 

''Errors are undetermined for this value. 
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Table 18. VV114 ISM Absorption Line Fit Data 



Ion 


Ao 


log{\fN/NH) 


Instrument 




Vc 


FWHM 




A 






A 






cm 


977.02 


-0.61 


lalif 


2.55±0.81 


5971±8 


350'' 








2blif 


2.39*" 


5990±5 


304±5 


oVI 


1031.93 


-1.13 


lalif 


1.52±0.10 


5684±25 


750'' 








2blif 


1.51±0.47 


5759±28 


783±196 


ClI 


1036.34 


-1.37 


lalif 


2.70±0.48 


5964±22 


746±134 








2blif 


2.54±0.17 


5962±17 


635±66 


ol 


988.73 


-1.56 


lalif 


1.59*' 


6084±312 


350'' 








2blif 


2.32±0.80 


6114±40 


448±37 


Nil 


1083.99 


-2.00 


Iblif 


2.25±0.09 


5975±10 


651±30 








2alif 


2.34±0.09 


5980±8 


645±26 


NlII 


989.80 


-2.03 


lalif 


2.33'' 


5873±11 


400±129 








2blif 


2.32±0.57 


5900±2 


395±37 


Ol 


1039.23 


-2.29 


lalif 


0.61±0.07 


6088±22 


310'' 








2blif 


0.53±0.10 


6110±23 


310±61 


Nl 


1134.41 


-2.57 


2alif 


0.35±0.06 


6021±27 


300±51 


Fell 


1096.88 


-2.96 


2alif 


0.15±0.05 


6055±33 


220'' 


SlII 


1012.50 


-3.24 


lalif 


0.30±0.09 


5964±24 


201±64 








2blif 


0.32±0.06 


6021±20 


211±40 


SilVs 


1128.33 




2alif 


0.13'' 


6024±31 


339±72 



^Using Vsys = 6016 kms ^ 

''Errors are undetermined for this value. 
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Table 19. O VI Data: Apparent Optical Depth Method 



Galaxy 


No VI 


b 


VQ VI — Vsya 




(cm-2) 


(kms^-'-) 


(kms^-'-) 


Haro 11 


8.0 X IQi* 


204 


-217 


VV 114 


2.2 X lO^s 


438 


-243 


SBS 0.335-0.52 


2.6 X 10" 


131 


-166 


NGC 1705 


2.1 X IQi* 


76 


-78 


NGC 1741 


6.3 X 10" 


228 


-127 


NGC 3310 


1.9 X 10^5 


232 


-325 


Haro 3 


7.5 X 10" 


149 


-102 


NGC 3690 


8.5 X 10" 


269 


-304 


NGC 4214 


4.0 X lO^* 


98 


-87 


Mrk 54 


1.3 X lO^s 


188 


-145 


NGC 7673 


1.9 X 10^5 


251 


-88 


NGC 7714 


1.7 X lO^s 


221 


-170 



Table 20. Line Data: Non-parametric Measurement 



c III c II N II 



Galaxy 


Weill 


vein 


Wen 


veil 


Wnii 


VNII 




(A) 


(kms-i) 


(A) 


(kms-l) 


(A) 


(kms-1) 


Haro 11 


1.74 


-253 


1.16 


-118 


0.92 


-113 


VV114 


3.37 


-212 


2.51 


-94 


2.11 


-84 


NGC 1140 


1.28 


-62 


0.87 


-45 


0.66 


6 


SBS 0335-052 






0.47 


-9 


0.28 


34 


Tol 0440-381 


1.28 


-151 


1.15 


-120 


0.71 


-82 


NGC 1705 


1.14 


-19 


0.87 


-38 


0.56 


-6 


NGC 1741 






1.27 


-41 


0.94 


-35 


IZW18 


0.61 


-137 










NGC 3310 


2.05 


-254 


2.54 


-186 


1.86 


-93 


Haro 3 


0.97 


-27 


0.74 


-49 


0.74 


49 


NGC3690 


1.91 


-153 


2.14 


-67 


1.80 


-77 


NGC 4214 


1.05 


-71 


0.93 


-18 


0.60 


-13 


Mrk 54 


2.15 


-167 


2.65 


11 


1.42 


-18 


IRAS 19245-4140 


1.71 


-98 


0.82 


-77 


0.20 


-21 


NGC 7673 


2.36 


-129 


1.44 


-96 


1.29 


-20 


NGC 7714 


1.91 


-46 


2.09 


-24 


2.69 


-67 



Note. — 



These measurements were made using the IRAF splot command "e" . 
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Table 21. Correlations of Outflow Speeds with Galaxy Properties 



Galaxy Parameter Outflow Measure Tq Pnull 



LiR - 






-vovi 


-0.636 


0.00644 


Lk 




Vsys 


- VOVI 


-0.564 


0.0158 






'^sys 


- vein 


-0.516 


0.0101 


Lk 




Vsys 


- Vein 


-0.495 


0.0138 



Note. — Col 3. measures the strength of the correlation 
between the parameters in the first two columns, using the 
value Ta from the Kendall rank test. Col. 4 specifies the 
probability of the null hypothesis that the two parameters are 
uncorrelated. 



Table 22. Correlations of Equivalent Widths with Galaxy Properties 



Galaxy Parameter 


Line 


Ta 


Pnull 


log{0/H) + 


12 


EQWo/ 


0.748 


0.000710 


Lk 




EQWo/ 


0.689 


0.00556 


LiR + Ljjy 




EQWoj 


0.600 


0.0157 


log{0/H) + 


12 


EQWc// 


0.694 


0.000960 


Lk 




EQWcj/ 


0.718 


0.000634 


LiR + LjjY 




EQWc// 


0.564 


0.00727 


log(0/H) + 


12 


EQWs/// 


0.442 


0.0753 


Lk 




EQWs/// 


0.697 


0.00283 


LiR + Luv 




EQWs/// 


0.624 


0.00756 


log(0/H) + 


12 


EQWc/// 


0.379 


0.0712 


Lk 




EQWc/// 


0.615 


0.00341 


LiR + Luv 




EQWc/// 


0.590 


0.00510 



Note. — Col 3. measures the strength of the corre- 
lation between the parameters in the first two columns, 

using the value Ta from the Kendall rank test. Col. 4 
specifies the probability of the null hypothesis that the 
two parameters are uncorrelated. 
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Table 23. Lyman Limit Escape Fraction 



Galaxy 


MW E(B-V) 


F 

1150A 


F 

1150A,corr 


F . 

900A 


F = 

900A,corr 


^iisoA+^ra 


L = 

900A 


fesc,abs 


fcsc,rcl 




(mag) 


(erg cm~^ 






(erg cm-^ 




L© 


L© 


(%) 


(%) 


Haroll 


0.011 


5.2 X 10"" 


5.9 X 10- 


-14 


< 1.3 X 10-1^ 


< 1.6 X 10-15 


2.1 X 10" 


4.6 X 10* 


<1.6 


<ii 


Mrk54 


0.015 


3.2 X 10-" 


3.8 X 10- 


-14 


< 6.5 X 10-1** 


< 8.9 X 10-1*5 


1.7 X 10" 


1.2 X 10** 


<5.0 


<9 


NGC 1741 


0.051 


5.2 X 10-" 


9.2 X 10" 


-14 


< 1.3 X 10-^'' 


< 3.7 X 10-15 


4.9 X IQio 


4.4 X 10** 


<6.3 


<16 


NGC 3690 


0.017 


7.4 X 10-" 


8.9 X 10- 


-14 


< 1.6 X IQ-^'' 


< 2.3 X 10-15 


6.7 X 10" 


1.6 X 10® 


<0.2 


<10 


VV114 


0.016 


7.6 X 10-" 


9.1 X 10- 


-14 


< 4.2 X 10-1^ 


< 5.8 X 10-15 


5.4 X 10" 


1.6 X 10'^ 


<2.0 


<26 



Note. — We have estimated an upper limit to the flux below the Lyman limit for the galaxies with S/N>5 in the LiF lA channel. See 
Section l3.5.1l for a discussion of the derived values in this table. We find fesc < 7% for all five galaxies. 



Table 24. Lyman Limit: C II Limits on the Escape Fraction 







P.F. 




Slab 




Galaxy 


fcsc.C 11 


fesc,abs 










(%) 


(%) 


(Z©) 


(kms-i) 




Haroll 


<8 


<0.6 


0.17 


225 


281.2 


Haro 3 


< 0.5 


< 0.1 


0.46 


141 


oo 


IRAS 19245-4140 


<21 


<11 


0.22 


163 


100.0 


IZwl8 


<10 




0.03 


111 


633.9 


NGC 1140 


<4 


<1.4 


0.22 


181 


229.0 


NGC 1705 


<1 


<0.8 


0.22 


155 


280.6 


NGC 1741 


<1 


<0.2 


0.22 


161 


291.5 


NGC 3310 


<5 


<0.5 


2.19 


172 


20.3 


NGC 3690 


< 1 


< 0.1 


1.38 


477 


oo 


NGC 7673 


< 10 


< 1.7 


0.69 


229 


oo 


NGC 7714 


<3 


<0.1 


1.10 


323 


88.8 


BBS 0335-052 


<18 




0.05 


66 


208.1 


Tol 0440-381 


<4 


<1.6 


0.34 


211 


172.4 


VV114 


<2 


<0.1 


0.87 


419 


161.9 


NGC 4214 


<1 


<0.2 


0.35 


66 


75.4 



Note. — We have conservatively estimated the continuum level and 
found the lowest convincing escape fraction in C II at line center. We 
then estimate fcsc.abs using the picket fence (P.F.) model for escaping 
radiation. We also find extremely high optical depths (> 20) below the 
Lyman limit assuming the slab model. 

''Derived using the gas p hase oxygen abundanc es in Table[l]and a solar 
C/O abundance ratio from lXsplund et al.l ll2005l) . 

''From gaussian fits in Tables 131-1181 
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Wavelength Shift [A] 



Fig. 1. — Histogram of all of the individual wavelength offsets used to correct for zero point 
shifts in the wavelength cahbrations. The vast majority of the corrections are less than 0.02 
A (6 km s~^). The overlaid arrows are 5"in length. 
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Fig. 2. — UV or FUV images with FUSE apeitme overlay for Haroll, VV114, 
NGC1140, SBS 0335-052, NGC1705, NGC1741, IZwl8, NGC3310, and Haro3. The 
FUSE LWRS aperture has been overlaid at the nominal position and orientation for the 
Ft/^i? observations. The image of Haro3 is a GALEX NUV observation at significantly 
lower spatial resolution than observed in the other panels. The overlaid arrows are 5"in 
length. 




Fig. 3.— UV images with aperture overlay for NGC3690, NGC4214, Mrk54, 

IRAS 19245-4140, NGC7673, and NGC7714. The FUSE LWRS aperture has been over- 
laid at the nominal position and orientation for the F?7S'i? observations. The images of 
NGC4214 {GALEX NUV) and Mrk54 (SDSS U-band) are at significantly lower spatial 
resolutions than observed in the other panels. 



-57- 




Fig. 4. — Spectra of Haro 11 from the FUSE SiC 2A, LiF lA, and LiF 2A channels. Locations 
of prominent ISM, stellar photospheric (dashed lines), Milky Way, and airglow (®) features 
have been identified. 
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Fig. 5.— Spectra of W 114 from the FUSE SiC 2A, LiF lA, and LiF 2A channels. Locations 
of prominent ISM, stellar photospheric (dashed lines), Milky Way, and airglow (®) features 
have been identified. 
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Fig. 6.— Spectra of NGC1140 from the FUSE SiC 2A, LiF lA, and LiF 2A channels. 
Locations of prominent ISM, stellar photospheric (dashed lines), Milky Way, and airglow 
(©) features have been identified. 
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Fig. 7.— Spectra of SBS 0335-052 from the FUSE SiC 2A, LiF lA, and LiF 2A channels. 
Locations of prominent ISM, stellar photospheric (dashed lines), Milky Way, and airglow 
(©) features have been identified. 
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Fig. 8.— Spectra of Tol 0440-381 from the FUSE SiC 2A, LiF lA, and LiF 2A channels. 
Locations of prominent ISM, stellar photospheric (dashed lines), Milky Way, and airglow 
(©) features have been identified. 
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Fig. 9.— Spectra of NGC 1705 from the FUSE SiC 2A, LiF lA, and LiF 2A channels. 
Locations of prominent ISM, stellar photospheric (dashed lines), Milky Way, and airglow 
(©) features have been identified. 
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Fig. 10.— Spectra of NGC 1741 from the FUSE SiC 2A, LiF lA, and LiF 2A channels. 
Locations of prominent ISM, stellar photospheric (dashed lines), Milky Way, and airglow 
(©) features have been identified. 
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Fig. 11.— Spectra of IZW18 from the FUSE SiC 2A, LiF lA, and LiF 2A channels. Lo- 
cations of prominent ISM, stellar photospheric (dashed lines), Milky Way, and airglow (®) 
features have been identified. 
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Fig. 12.— Spectra of NGC3310 from the FUSE SiC 2A, LiF lA, and LiF 2A channels. 
Locations of prominent ISM, stellar photospheric (dashed lines), Milky Way, and airglow 
(©) features have been identified. 
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Fig. 13. — Spectra of Haro 3 from the FUSE SiC 2A, LiF lA, and LiF 2A channels. Locations 
of prominent ISM, stellar photospheric (dashed lines), Milky Way, and airglow (®) features 
have been identified. 
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Fig. 14.— Spectra of NGC 3690 from the FUSE SiC 2A, LiF lA, and LiF 2A channels. 
Locations of prominent ISM, stellar photospheric (dashed lines), Milky Way, and airglow 
(©) features have been identified. 
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Fig. 15.— Spectra of NGC4214 from the FUSE SiC 2A, LiF lA, and LiF 2A channels. 
Locations of prominent ISM, stellar photospheric (dashed lines), Milky Way, and airglow 
(©) features have been identified. 
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Fig. 16.— Spectra of Mrk 54 from the FUSE SiC 2A, LiF lA, and LiF 2A channels. Locations 
of prominent ISM, stellar photospheric (dashed lines), Milky Way, and airglow (®) features 
have been identified. 
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Fig. 17.— Spectra of IRAS 19245-4140 from the FUSE SiC 2A, LiF lA, and LiF 2A channels. 
Locations of prominent ISM, stellar photospheric (dashed lines), Milky Way, and airglow (®) 
features have been identified. 
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Fig. 18.— Spectra of NGC7673 from the FUSE SiC 2A, LiF lA, and LiF 2A channels. 
Locations of prominent ISM, stellar photospheric (dashed lines), Milky Way, and airglow 
(©) features have been identified. 
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Fig. 19.— Spectra of NGC7714 from the FUSE SiC 2A, LiF lA, and LiF 2A channels. 
Locations of prominent ISM, stellar photospheric (dashed lines), Milky Way, and airglow 
(©) features have been identified. 
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Fig. 20. — Prominent spectral absorption lines plotted for Haroll. Note: O IA989 is blended 
with N IIIA990. 
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Fig. 21. — Prominent spectral absorption lines plotted for VV114. Note: O IA989 is blended 
with N IIIA990. 
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Fig. 22. — Prominent spectral absorption lines plotted for NGC1140. Note: O IA989 is 
blended with N IIIA990. 
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Fig. 23. — Prominent spectral absorption lines plotted for SBS0335-052. Note: O IA989 is 
blended with N IIIA990. 
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Fig. 24. — Prominent spectral absorption lines plotted for Tol0440-381. Note: O IA989 is 
blended with N IIIA990. 
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Fig. 25. — Prominent spectral absorption lines plotted for NGC1705. Note: O IA989 is 
blended with N IIIA990. 
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Fig. 26. — Prominent spectral absorption lines plotted for NGC1741. Note: O IA989 is 
blended with N IIIA990. 
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Fig. 27. — Prominent spectral absorption lines plotted for IZW18. Note: O IA989 is blended 
with N IIIA990. 
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Fig. 28. — Prominent spectral absorption lines plotted for NGC3310. Note: O IA989 is 
blended with N IIIA990. 
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Fig. 29. — Prominent spectral absorption lines plotted for Haro3. Note: O IA989 is blended 
with N IIIA990. 
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Fig. 30. — Prominent spectral absorption lines plotted for NGC3690. Note: O IA989 is 
blended with N IIIA990. 
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Fig. 31. — Prominent spectral absorption lines plotted for NGC4214. Note: O IA989 is 
blended with N IIIA990. 
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Fig. 32. — Prominent spectral absorption lines plotted for Mrk54. Note: O IA989 is blended 
with N IIIA990. 
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Fig. 33. — Prominent spectral absorption lines plotted for IRAS19245-4140. Note: O IA989 
is blended with N IIIA990. 
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Fig. 34. — Prominent spectral absorption lines plotted for NGC7673. Note: O IA989 is 
blended with N IIIA990. 
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Fig. 35. — Prominent spectral absorption lines plotted for NGC7714. Note: O IA989 is 
blended with N IIIA990. 
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Fig. 36.— Overlay of O IA976 on C IIIA977 for NGC 1140 and NGC 3310. The O I feature 
is significantly weaker than the C III profile and does not significantly affect our measured 
parameters for the C IIIA977 line. 
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Histogram of Outflow Velocities in Order of Increasing Ionization Energy 
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Fig. 37. — Velocity Outflow Histogram for Prominent Absorption Lines: Outflows are 
observed for all of the absorption features except O IA1039. Observed outflow velocities are 
significantly higher for O VI and C III. 
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Fig. 38. — Outflow Velocity Comparison: The plot on the left shows the C IIA1036 versus 
N IIA1084 velocities while the right displays the C IIIA977 versus O VIA1032 velocities. 
C II and O VI are generally more blueshifted than the N II and C III features respectively. 
The lines represents Vx = Vy for both plots. NGC 3310 (which has unusually large outflow 
velocities for its mass and star formation rate) is denoted by a red circle. 



C II versus S III Outflow Velocities 




Fig. 39. — C IIA1036 versus S IIIA1012 velocities. While S III has a significantly larger 
ionization energy (23.3 eV) than C IIA1036 (11.3 eV), higher outflow velocities are observed 
in C II. The line represents vcn = vsni in both plots. NGC 3310 (which has unusually large 
outflow velocities for its mass and star formation rate) is denoted by a red circle. 
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Fig. 40. — Column density vs. line width for a wide variety of O VI absorption line systems 
including galactic disk and halo, high velocity clouds, starburst galaxies, and the IGM. The 
two dashed lines indicate the predictions for radiati vely cooling gas for ass umed temperatures 
of Tovi = 10^, 10^ K. This plot is originally from (iHeckman et al.l 120021 ). The star-forming 
galaxies from this work have been added and extend the relationship to much higher column 
densities. 
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Fig. 41. — The log of the O VI outflow velocities in km/s versus the log of the K-Band and 
IR+UV luminosities (in Lq). These are roughly proportional to the stellar mass and star 
formation rate respectively. NGC 3310 (which has unusually large outflow velocities for its 
mass and star formation rate) is denoted by a red circle. 
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Fig. 42. — Same as Figure 41, except that here we plot the 'maximum' outflow velocity. 
This is defined so that 90% (10%) of the absorbing column of gas is to the red (blue) of the 
velocity. See text for details. NGC 3310 (which has unusually large outflow velocities for its 
mass and star formation rate) is denoted by a red circle. 
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Fig. 43. — The log of the C III outflow velocities in km/s versus the log of the K-Band and 
IR+UV luminosities (in Lq). These are roughly proportional to the stellar mass and star 
formation rate respectively. NGC 3310 (which has unusually large outflow velocities for its 
mass and star formation rate) is denoted by a red circle. 
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Fig. 44. — O VI and C III outflow velocities versus (Ljr + Luv)/Lk. While these plots 
have a large degree of scatter, the higher SFR/Mass galaxies tend to have higher outflow 
velocities. NGC 3310 (which has unusually large outflow velocities for its mass and star 
formation rate) is denoted by a red circle. 
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Fig. 45. — Equivalent widths of O I, C II, S III, and C III compared to gas phase metaUicity, 
K-band luminosity, and IR+UV luminosity. For the lower ionic states (O I & C II) the gas 
phase abundance and stellar mass track the strength of the lines. The higher ionization 
states however poorly match the metaUicity and instead are correlated with the SFR and 
stellar mass. 
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Fig. 46. — This figure shows the Lyman continuum regions of five galaxies in order of 
increasing redshift. We have plotted only the night data in order to minimize contamination 
from airglow and earth limb emission. There is no convincing evidence of Lyman continuum 
emission in these spectra. The spectra are generally very noisy in this wavelength region 
and any possible signal is dominated by the background uncertainties and airglow emission. 
These plots do show earth limb and airglow contamination from H I, O I, N I, and N2 (for a 
better example see Figure Wl} . 
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Fig. 47. — This is identical to Figure HE] but for an additional four galaxies. No convincing 
evidence of Lyman continuum emission is observed in these galaxies. Earth limb and airglow 
emission however are present from H I, O I, N I, and N2. This is easily confirmed by compar- 
ing the specta, particularly VV 114 and Tol 0440-381 which both have strong contamination 
from airglow in the Lyman continuum regions. These features however are present in all 
eight spectra in the two figures. We have included both the night and day time data for 
VV114 due to the small percentage of night time data available for that target. 
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Fig. 48. — Plo ts of LWRS SiC IB nig ht airglow (top) and earth limb (bottom) spectra. Data 
originally from iFeldman et al.l (120011 ). Though the night airglow spectrum is very noisy the 
two spectra show the wealth of geo-coronal emission features that impact data analysis in 
these wavelength regions. 



